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ELECTRICAL NOTES 
By F. B. PIDDUCK (Oxford) 
1X. OSCILLATIONS WITH A SPLIT-ANODE MAGNETRON 
[Received 27 June 1936; in revised form 19 August 1936] 


1. Introduction 

In applying the methods of my last notet to split-anode magnetrons 
I have continued to neglect the inequality of emission velocity and 
to assume that most of the electrons are caught on the filament 
when they return. This restricts us to emissions of the order of a 
microampere, and also cuts across another assumption, necessary 
for a simple theory, that the paths of the electrons are nearly circles. 
I therefore regard this note as answering a qualitative question, 
whether known electrical principles can account for the maintenance 
of any short-waved oscillations by a split-anode magnetron, by taking 
a strict theory and stretching the conditions of its application. The 
result appears, within its limits, to be encouraging. 


2. Motion of electrons between a filament and a split anode 

The currents are supposed to be so small that the electric and 
magnetic field of the electrons themselves has no appreciable effect 
on their motion, the electrons being taken to move in a field of scalar 
electromagnetic potential V, on which a uniform magnetic field H is 
superposed parallel to the filament. Let —e be the charge and m the 
mass of an electron, 


w = He/me. (1) 


The equations of motion are 


. wie 
L+wy = ——, 
m Cx 


Let the motion in an assigned static field be known, and let x(t), y(é) 
be the coordinates at time ¢ of an electron leaving a given point at 
time ¢, with given velocities in the plane of x,y. Let V now denote 
the scalar electromagnetic potential of a small field superposed on 
the static field, so that V is a given function of x, y and ¢, and let 
X(t), Y(t) be the result of substituting x(t), y(t) for x, y in éV/éx, 


+ See above, pp. 210-213. 
R 
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éV /éy respectively. Then if x, y are now the small changes of 
coordinate at time f, 


+ ij—iw(é+iy) = (e/m)| X(t)+7Y (6)] 
and x = y= = y = 0 whent = ¢,. Hence 


t 
xtiy = am | [X(r)-+-¢¥ (7) |[1—exp i(wt—wr) | dr. (2) 
mw : 

Let a be the radius of the filament, small and sometimes taken 
as zero, b the radius of the anode, and let the anode be divided by 
a fine cut at polar angles +47. Let all the electrons be emitted 
normally with the same velocity. Then in the steady motion, with 
the same constant potential applied to both halves of the anode, all 
are turned back at the same radius, and we shall examine the effects 
produced when that radius is nearly 6. 

Let the disturbance consist of potentials +}V applied to the right- 
and left-hand halves of the anode, where V is in general a small 
function of ¢, but such that the wave-length in air of the correspond- 
ing electric oscillations is much greater than b. The effect of the 
filament on the field so added will be neglected, so that the added 


potential at a point z = x+y inside the circle |z| = 6 is the real 
part of the function 
: b+iz 
W = — Tog sand 
oh —iz 
Thus X(t)+iY (t) = 26V /n[b?+-2*(t)?]. (3) 


For an electron leaving the filament at an angle 0, at time fo, 


x 


(t) = tb exp 10,[ 1— exp t(wt—wty) |. 
Let V = Acosvt. (4) 


Then from (2), the electron is at time ¢ at the point 


z= 3tb&(1—n)+bgAET(E, 0), (5) 
where 
€ = exp 10,, ” = exp i(wt—wty), g = 8e/mmw*b? (6) 
wt — wt, 





Sn ena 
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Writing x = exp(—i¢), cos(vtj+vd/w) is a function of 2, which we 
call F(x), and 1/n 

a EF (x)(—nx+1) dx 
Ke) = | oe wee ”" 
1 


taken along an arc of the circle |z| = 1 below the real axis. 

At a point on the boundary of the cloud of electrons, x and y are 
unaltered by certain small changes of 6, and t). Hence the condition 
to be satisfied on the boundary is @(x, y)/@(A,t)) = 0, or 

O(z, z " “a 
ag, n) 
in which €*, n* are to be replaced by €-1, n-1. It is easy to see that 
wt—wt, = 7 in the steady motion, so that 1+ 7 is small and imaginary. 
A slight error in 7 does not affect the radial distance of an electron 
on the boundary, and it follows from (5) that that distance is given by 


(8) 





r/b = 1+gA Ximaginary part of J(€,—1). (9) 
To calculate 1+ 7, write IJ for I(€,) and IJ-+ for I(€-4,y-1). Then 
J*—J~-* is an integral round the whole circle |x| = 1, so that z* is 


expressed easily in terms of € and 4. The integrand in (7) is rational 
if v/w is an integer, but its residue at the pole x = 0 is different in 
each case. More generally, J is the integral of a rational function if 
v/w is commensurable. 


3. Boundary of the cloud of electrons when v = 0 
Here F(x) = 1, J* = J-1+ 2imé(é*—4)-1 and (8) becomes, correct 
to terms of the first order of small quantities, 
1+» = igA[I+€0l/0E+20I/0n+I7+€0I 1/04 +- 

+ 20I-1/an-1+ 4in€3(€2@—4)—], 
where 7 is put equal to —1 in the small terms after differentiation. 
The conjugate complex of J+£é0J/éé+2éI/én differs from the corre- 
sponding —1 ere by 4imé3(é°—4)-*. Write 


&(x+1) dx 
‘* fia a—1—2€)(a—1+-2€)’ 
> (10) 


(w—1—2£)%(a—1+26)*  * 
Then 1+ = 2igA xreal part of (—J—K). 





x= [> aa (w+-3)(~—1)?+-4€*(~—1)] da 
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If we put J = J,+7J,,..., where J, and J, are real functions of €, the 
polar coordinates of the electron on the boundary which left .the 
filament at angle @, are given by 
r/b = 1+gAJ,, 6 =437+0)+gAK,. (11) 
We find 
4(17—8 cos 20,)J, = (—10sin &)+4 sin 20,)log(2 cos $65)— 
—(10sin #,+ 4 sin 20,)log(+ 2 sin $6))— 
— 6 cos 6,(0,— 27) + (—9-+3 cos 6)+-6 cos 26)7, 
4(17—8 cos 20,)*K, 


4 cos §)?(— 10—18 cos 0,—8 cos 20,)log(2 cos $6,)+- 





== (5 
+ (5+4 cos ,)?(10—18 cos 0,+-8 cos 26,)log(+ 2 sin $0)+ 
+ 24(17—8 cos 26,)cos 0,+ (246 sin 0,— 104 sin 30))(@,—2a)+ 
+ (5+ 4 cos 6,)?(—7 sin 0)+-4 sin 26,)z, 
where —z < 6, < wand the upper or lower signs are taken according 
as @, is positive or negative. The functions are tabulated below from 
0 to 180°, after which they repeat with reversed sign. 


8, J. Ky 8, J. Ky 8, J; Ky 

0 1-047 0-590 | 60° | —0-240 —0-759 | 120 0-689 | —0-758 
10 1-054 0-893 70° | —0°369 | —0-856 | 130° | —0-727 | —0-718 
20 0-845 0-744 80° —0-466 0-880 | 140° —0-766 —0-666 
30 0-457 0-269 90° | —0-541 0-880 150 —0O-798 | —0-622 
40 0-161  —0-214 |) 100° —0-600 -O-851 | 160° | —0°844 —0-586 
50 — 0-068 -0-556 | 110° | —0-648 0-820 | 170 -0-903 —0-568 


It is well known that a magnetron can maintain slow oscillations 
explained by the characteristic falling over part of its range. The 
table illustrates this to some extent. Let A be positive, so that the 
right-hand sector is at the higher potential. Since J, is greatest when 
6, is about 6°, the radial extension of the cloud is greatest near 
6 = 96°, so that an anode slightly larger than the undisturbed cloud 
is grazed first at a point on the negative sector. If the anode is the 
exact size of the undisturbed cloud, an infinitely small potential 
brings electrons between 6, = 0 and 47° to the negative sector and 
electrons between 227° and 360° to the positive sector, that is more 
to the positive sector. But since K, = 0-59 when 6) = 0, the 
electrons falling on the negative sector range from 0, = —0-59gA 
to 47° when A is small but finite, so that the range increases with A 
and the characteristic falls just after its discontinuous rise. 
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‘4. Boundary of the cloud when v = w 
Here I(€, —1) = —}3[K exp twt,+ Lexp(—iwt,)], where 
-1 
—E&(x+1) dx —&(x+1) dx 
" “Li wane &— | eee 
The same letter K has been used since the integrals in (10) and 
(12) are equal. Hence from (9) the radial distance of an electron on 
the boundary is given by 
r/b = 1+gApcos(wt,—a), (13) 
where pcosa = —}(K,+L,), psina = —}(K,—JL,). We find 
4(17—8 cos 26))?K, 
= (5—4 cos 4,)?(14 sin #,+ 8 sin 26,)log(2 cos $6))+ 
+(5+-4 cos 6,)?(14 sin #,—8 sin 26,)log(+-2 sin $0)— 
—40(17—8 cos 26,)sin 6,+-(106 cos 8,— 88 cos 36)(0g— 2m) + 
+(5+4 cos 6))?(5—9 cos @,+ 4 cos 26,)7, 
4L, = —2(1+ cos 6))log(2 cos $6)+- 
+ 2(1—cos 6))log(+ 2 sin $6) — 2 sin 6) -+-sin 4% 7, 
4L, = —2sin @ log(2 cos $4))— 
— 2sin 0, log(+ 2 sin $0,)-+ 2 cos 4 @,+(1—cos 0))z. 
After 180°, p repeats and a is increased by 180°. 





0 0-664 ; 60° | 0-378 49° 24’ 120° | 0-541 | 119° 30’ 
10° 0-766 j 70° | 0-385 72°18’ | 130°| 0-565 121° 6’ 
20 0-746 35° 6 80° | 0-410 90° 48’ 140° | 0-581 | 122° 
30 0-630 90° | 0-448 | 100° 42’ 150° | 0-599 | 121° 48’ 
40 0-503 j | 100° | 0-481 | 108° 54’ | 160° 0-617 120° 24’ 
50° | 0-410 110° | 0-517 | 114° 18’ | 170° 0-641 | 116° 48’ 


8 p 0, p o —% ; p a 


If the radius of the anode is exactly b, electrons are caught on it 
when the point 65, ¢) lies between the curves wi,y—a = +42. The 
range of 0, is 180° for most values of t), but falls short of it when 
15° < wt) < 28° and 195° < wt, < 208°, on account of the strong 
field near the edges of the sectors. It is simpler to ignore the defect 
and take the range to be 180° throughout. Then half of the electrons 
emitted between times ¢) and t)+dt, are caught on the anode, which 
is equivalent to supposing that a thin semicircular sector of positive 
charge }cJI dt, is emitted by the anode at time ¢,+7/w and moves 
towards the filament so that it remains nearly semicircular, the ends 
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moving on ares of the circle z = }ib€(1—7). The initial semicircle 
moves round the anode as f, increases, 6, lying for example between 
30° 42’ and 210° 42’ when wt, = 60° and between 53° and 233° when 
wt, = 120°. The aggregate of all the contracting semicircles gives 
rise to a current as explained in § 5. 


5. Calculation of current 

Let e be a set of charges inside the magnetron at time ¢, and let 
S,, S, be the two sectors of the anode, S, extending from polar angles 
47 to 37 and S, from —47 to 3x. The filament S, is supposed to be 
too small to affect the charges e,, e, induced by e in S,, 8, when at 
some common potential, or the mutual influence of S,, S,. Then 
€: +e. = —e, Gor = Gog ANA Gy, = Yoo. The sum of the charges on S, 
and S, is —e—2q9, Y+(d41+912)(4+) and the difference is 


a 


€y—€2+ (F11— Mi2)(h—hh). 

If %) and V,+J, are constant, the sum of the charges is constant, and 

the rate of increase of the charge on S, can be regarded as a current 
from S, to S, of amount 

KdV | de, 


14 
dt ‘ cdt’ ome 


where V is the potential of S, above S, and K the capacity of a con- 


denser whose plates are S,, S5. 

For a charge e at a complex point z, the image is a charge —e at 
the inverse point b?z*-1 and the complex potential at a point b exp 70 
on the anode is 

V+iU = —2elog(bexp 10—z)-+ 2e log(b exp 10—b?z*-), 
where e is reckoned temporarily as a charge per cm. Since —4ize, is 
the difference between the values of V+-iU when @ = $a and —}$z, 
2ime,/e = —log(th—z)+log(—1b—z)+ 
+log(ib—b*z*-1)—log(—ib—b?z*-1), 

It is sufficient to take z from the steady motion, thatisz = 3ib€(1—»). 
Since dyn/dt = iwy, 

2rde, 


wedt 


ae 
2+€—fn 2—§+&4 

ees. a Le 

(n—1—2€-1)(n—1+-2-*) 


4 2én ; __ 2&m 


(1—n)(1— + 2&m) © (1—n)(1——2€n) 
4€y-1 
+ (qA—1— 24)(q— 1428)" 
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For a charge e distributed over a semicircle and moving inwards 
as described in § 4, replace e by the charge ed6,/7 between 0, and 
6,+-d6), writing it as iedé-1/7é-" in the first term and —iedé/7é in 
the second and integrating. Hence the current at time ¢ due to such 
a semicircle, leaving the filament between polar angles @, and 6,+-7 
at time t), is RwI dt,/7*, where 

9-1-5 

saat ee) 

and the correct branch of the logarithm must be taken in each case 
to represent the continuous passage of @, from one end of its range to 
the other. The function R is tabulated on p. 248 for 10° intervals of wt 
and wt—wt,, where wi—wt, goes from 180° to 360° and R repeats with 
reversed sign when w¢ is increased by 180°. The columns of the table 
have the same relation to the total current that the ordinates of a 
curve have to its area, and are treated by Simpson’s rule so as to 
give that current. The result is shown in the following table: 


R = real part of -_ los( — 
n—1 








Convection | Convection Convection 

wt current/I wt |  current/I wt current/I | 

0 —0-222 | 60° — 0-038 120° 0-198 
10° — 0-200 70° 0-015 130° | 0-225 
20° —0-176 80° 0-049 140° 0-242 
30° —0-139 | 90° 0-096 ~—-:150° 0-248 | 
40° —O-117 100° 0-129 160° 0-248 
50° —0-068 | 110° 0-169 170° 0-239 | 








Write 4(wt) for the function just tabulated in the range 0 < wt < z. 
Then if it is expanded in a Fourier’s series > (a,, cos nwt-+-b, sin nwt) 
in the range 0 < wt < 27, its fundamental is the real part of 


(a,—1ib, exp twt, where 
9 7T 
ih ee | $(0)exp(—i8) dd. 
7 
0 


We find by Simpson’s rule that the fundamental term in the con- 
vection current is the real part of 
(—0-233—0-098i)/ exp iwt. (15) 


6. Oscillations on Lecher wires 

We use the theory (reasonably valid when the frequency is high and 
the resistance low) which works with self-inductance L, resistance 
R and capacity C per cm., where R is the resistance per cm. of either 
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wire at the frequency v of the oscillations, and LCc? = 1. Let the 
wires be joined at 2 = 0 to the sectors of the anode and bridged 
across at 2 = /, and let the high potential be applied to the anode 
through large chokes. Then there is a solution in which the current 
at distance x from the magnetron is 

[P exp itka+ Q exp(—ikzx) lexp it 
and the potential 
(k/vC)[P exp tka—Q exp(—ikx)]exp ut, 

where k? = v?/c?—2iv RC; or if R is small 

k = v/e—icRC. (16) 
The potential vanishes when 2 = /, and hence the current at « = 0, 
that is the current through the magnetron, is 

(ivCV /k)cot kl, (17) 
where the signs are chosen so that the current flows from S, to S, and 
the real parts of the complex quantities are taken in the physical 
interpretation. 

Let the magnetron be working as described in § 4 and 5, that is 
let v = w and let the radius of the anode be b. Then we have from 
(4), (14), (15), (16) and (17) 

| ivk A + (—0-233—0-0982) J |(vp—tc? RC) — 
—icvC A (cot vl/e+-ic RCI cosec*vl/c) = 0. 
Neglecting products of R and J and equating real and imaginary 
parts to zero, 
A(vK—cC cot vl/c)—0-098] = 0, 
c?RCA(K-+ Cl cosec*vl/c)—0-233] = 0. 


The zero approximation R = J = 0 gives the natural frequencies 
vtanvl/ce = cC/K. A next approximation gives frequency v+6r, 
where dv = 0-42R/L. The amplitude of potential maintained is 


A = 0-2331L/R(K+Clcosec*rl/c). (18) 


The condition of maintenance, that (18) makes A positive, is satisfied. 

Take for example a magnetron with anode 4 cm. in diameter, 
magnetic field 50 gausses, emission 3x 10-* amperes, attached to 
Lecher wires of copper, each 1 mm. in diameter and 2 cm. apart. 
The capacity of the magnetron itself is taken as 1 cm. Then 
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v = w = 8-9 108, the resistance of each wire at that frequency is 
0-009 ohms per cm., the steady anode potential 225 volts, 1 = 39 em. 
for the lowest mode, A = 0-246 volts and dv/v = 1/3500. 


7. Effect of the filament on returning electrons 
We shall discuss the case v = w, for which 
1/y 
expiwt, [ &(—na+1) dx 
2 x?(x— 1— 2£)(~a— 1+-2€) 


© \ 
1 


I(é, n) ‘ 


1/7 

_ exp(- tut») " &(—nx-+1) dx 

2 J) (w—1—2€)(~—1+2€) 
1 


If wi—wt, = 27, integration is round the whole circle. Thus the 
electron would be at the point 
z —imbgA€*(1-+-4€*)(1— 4€*) "exp iwt, 
at time ¢,-+27/w, if it were not for the repulsion of the filament, and 
would have a complex velocity w nearly equal to }bwé. This electron 
approaches the filament with a moment of velocity 
h 1i(zw*—z*w) 
kab*gwA real part of €(1-+4€7)(1—4£?)~*exp iw, 
round it. The greatest numerical value of h, considered as a function 
of ty, is 


Lirb?qwA |(1+-4€?)/(1—4€?)?| = 8eA(17+-8 cos 26,)'/maw(17—8 cos 26), 


so that the greatest value of h for any returning electron is 20eA /9mw. 
If A = 0, finite velocity of emission will make some of the electrons 
(those, at least, which leave the filament normally) penetrate it and 
be caught on return. It remains to find whether, with any given 
amplitude of oscillating potential, most of the electrons are caught or 
not. The theory is simplified by assuming normal emission with a 
fixed velocity, corresponding to potential B. 

Let an electron at a large distance from the filament be moving 
directly towards it with such a velocity that its kinetic energy on 
arriving at the filament 7 = a is eB, and let it be displaced laterally 
so that the moment of its velocity becomes h. Then 

im (#2-+-1262)— el, log r/a ~ log b/a = const., 


where 726 = h and the constant is the same as for the undisplaced 
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electron, that is eB. If the displaced electron just reaches the fila- 
ment before turning, mh?/2a2 = eB. Since the high potential VY, is 
approximately equal to mwb?/8e, the least radius of filament which 


will catch all the returning electrons is given by 


8la*V, B = 25b7A?. (19) 
If ) = 225, A = 0-246 and B = 0-1, a/b = 1/38. This ratio is too 
large for the paths to be accurate circles, so that the theory is only 
qualitative at the best. 








YXPANSIONS BY THE FRACTIONAL CALCULUS 
By W. FABIAN (London) 


[Received 4 April 1936] 

1. Introduction 
d 
lx 
d A ] ae f™(a) 
x S (-)- ghte, l 
(,] J(@) (A) PM n! (A n) (1) 


A 
f(x) for real values of A by the series 


BERNOULLI* defines ( 
f 


Bourlet+ and Davis} have shown that this series can be obtained 


from Riemann’s definition 


d\- 1 fd\f 
= ; (x t)A ¥-lf(t) dt. 
(7) I(x) D(A = | I 
0 
where y is the least non-negative integer such that A+y > 0. With 
Riemann’s definition, which I have adopted in my papers on the 


fd\-A/d\., ; 
BS (--} J(%) I(x). 


This suggests that an expansion of f(a#) can be obtained under appro- 
, are : : ah 2 ee 
priate conditions by applying Bernoulli's formula to {|{—} f(x). Suit- 

ax 


Fractional Caleulus,* 


able conditions, under which a function can be expanded in this 
manner, will be determined in this paper. 
2. Fractional Integration 
A Ath integral of f(z) along a simple curve / is defined by 
Dr f | 
: (z—t)+7-f(t) dt, 
PiA+y). 


a 


D-X(1,,) f(z) 


* Pincherle, Mem. d. R. Acad. di Bologna (Sci.) (5) 9 (1902), 745, gives 
Bernoulli’s form of the fractional differential coefficient. 

+ Bourlet, Ann. de U’ Ecole Normal (3) 14 (1897), 154. 

t Davis, American J. of Math. 46 (1924), 95-109. 

Riemann, Ges. Werke (1876), xix. 331-44. 

** Fabian, Phil. Mag. (7) 20 (1935), 781-9; 21 (1936), 274-80. American 
J. of Math. and Phys. 15 (1936). Math. Gaz. 20 (1936), 88-92 ; 249-253. 

Expansion (1) is established here for the complex variable and complex 


values of A. 
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where y is the least non-negative integer such that R(A)+y > 0; 
D stands for d/dz, and the integration and differentiation are along 1. 

I have shown elsewhere* that, if f(z) is analytic throughout a 
bounded simply-connected region EL, which contains /, then each 
branch of D-(1,) f(z) is analytic inside EZ, except possibly at a. Under 
these conditions the value of each branch of D-(l,) f(z), given a and 
z, is independent of /. 


3. Expansions 
THEOREM 1. Let f(z) be analytic within a circle of centre z, which 
contains | in its interior. Let X be such that R(A) < 0, and let 


f(a) = f'(a) =... =f’a)=0 if RAY<- 


Then fo) ae 7 ake )”- J Trae —a)"-, 


provided* that the Taylor series for f\(z) at z converges on l uniformly 
to fy(z) 

Proof. By the properties of f\(z) given in § 2, the circle of con- 
vergence of the Taylor series for f\(z) at z, extends at least as far 
as a. Expanding therefore f\(t) at t = z) in a Taylor series, we have, 
since, by hypothesis, this series converges to f)(¢) uniformly on J up 
to a, 


A \(z i eins { —t)—- { Sr-u(%0) —-z n\ 
Dla) fx(@0) = r'(—A) | ¢ %—t) {> n! (t— Zp) ;* 


n=0 


—1)n Sr-nl@o) _(, 


—a)"-A, 


sai i ror} n! (n—A) 


By a previous theorem,~ 
D\(L,)fx(20) = f(%)- 

The proof of the theorem is completed. 

The expansion just obtained can also be proved to hold under 
different hypotheses. These are given in Theorem 3, for which the 

* Math. Gaz. 20 (1936), 249. 

+ fy_,(%) stands for D" AL) f (Zo)s and f,(z) stands for D-X(1,)f (2). This nota- 
tion will be used, when no ambiguity can arise. 

t Fabian, Phil. Mag. (7) 20 (1935), 781. 

The conditions f(a) = f’(a) = ... =f-(a) = 0 when R(A) < —1 are re- 
quired for this step of the proof. 
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following preliminary theorem is required: 


TuroreM 2. Let f(z) be analytic on l. Let X be any number, and let 


fle) =f(e)=..=fr*@)=0 Ff BA) < —1. 


Then, if 5 is any non-negative integer such that R(—A)+-6 > 0, 


n—1 
“i —l)m fy ml is 2- -a)™ at. 
m!(m - 


~ me. m! '(m—A-+-n)(z—a)"+8-A , 
Lad Z s!(m s)! D(m iw 1)P(n l ith a — ~ 


Z t ‘ 
eo? oma Pi soya dt 
for all positive integers n; all integrations are along 1, and the sum 


1 ™m 
S$ > is omitted if 6 0. 


_ 
m=+0s8s 0 


Proof. f(z) = D°D~(1,) fx(z), 
by a previous theorem,* 


Dp ¢ la 
: + t) A 1! )dr 
1'(d a 4 Sar) te 


n—I § 
Fey J eo tay DS Beamer | finale in| 


\m 


a 


t 
(on integrating { f\(r) dr by parts n times) 


Ds Dy m(2) (z—a)” oO A 
m! (m+5—A) 


De 


B+ t 
| —_tey © | 
+r »] (z—t) 4 frn(t) at ja 


which comes down to the form enunciated in the theorem, when the 


differentiations of the first » terms of this expansion are performed. 


* Fabian, Phil. Mag. (7) 20 (1935), 781. 





EXPANSIONS BY THE FRACTIONAL CALCULUS 


From the theorem just proved we can immediately deduce the 
following theorem, which gives the expansion in Theorem 1: 


THEOREM 3. Let f(z) be analytic on 1. Let X be any number, and let 
f@=f@=..=f7°@=0 ff RA)<-1. 


Then, for those values of z for which 


5 z t 
Poem | (env | tana) a 


tends to zero, and (z—a)"~fy_,,(z) remains bounded when the positive 
integer n tends to infinity, 
= fr—n(2) 
0 —])jr 24-0" _(z—a nd, 
fle) rw 2, Yee") 
all integrations being along 1, and 8 being a non-negative integer such 
that R(—A)+8 > 0. 





CONJUGATE MATRICES 
By A. R. RICHARDSON (Swansea) 
[Received 29 May 1936] 


1. CoNJUGATE matrices* are usually defined as a set satisfying the 
conditions: 

(i) they are commutative; 

(ii) their symmetric functions are scalar matrices which are equal 
to the corresponding symmetric functions of the scalar roots of the 
characteristic equation of one of the matrices; 

(iii) all the matrices have the same characteristic equation; 

(iv) all are of a specified order n. . 

P. Franklin} has generalized this by omitting condition (iii), whilst 
Sokolnikoff{ replaces (iii) by the condition that the matrices have 
the same reduced equation. Condition (iv) is not mentioned as such 
but is always implied. In some respects this theory resembles that 
of algebraic numbers but the analogy is very incomplete.§ For 
example, condition (iv) requires that the matrices involve the latent 
roots of the characteristic equation and also the nth roots of unity. 
In this note we show that, if in place of condition (iv) we require the 
matrices to be of the same order as the Galois group of the reduced 


equation, then there is a complete simple isomorphism between the 


fields of the original equation and that of its matrix roots. Also, if the 
elements of the matrix roots are to belong to the field of the coeffi- 
cients of the equation, then the order of the complete matrix-repre- 
sentation is the smallest possible. If the roots are represented by 
matrices of order higher than that of the Galois group, then the realm 
formed by adjoining them to the rational field may contain divisors 
of zero. Purely rational operations suffice to enable us to calculate 


the conjugate roots and the automorphisms of the Galois field. 


2. If the order of the matrix solutions of an equation is restricted 


to be that of the matrix coefficients, it appears that certain equations 


¥ Pp. yA 5. All references, unless otherwise stated, are to The Theory of 
Matrices, C. C. MaeDuffee (J. Springer, Berlin 1935). 
+ P. Franklin, pp. 24-5. 
I. S. Sokolnikoff, American J. of Math. 55 (1933), 167. 


§ A. A. Bennett, p. 24. 
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will have no solution, others a restricted number. For example, 


P. Franklin® shows that ° = 


0, 0 


ge 


(1) 


has no solution in matrices of order 2. Also it may be verified that 
2 
= 16 § 
00% 





has only one independent rational matrix solution of order 3. 
Normality is restored if the order of the matrix solution is raised. 
Thus (1) has the two rational conjugate solutions 
000 — 0 0O0T 
0 0 0 0 000 
100 ; —l1 000 
0: 0 0 —1 0 0. 
whilst (2) has the three rational conjugate solutions 
2000 | O QO —- — 
0000 
0000 
000 2 
010 
001 














0 
0 —1 0 —1 
0 o-1 O-—1 | 








Evidently the number of solutions is never reduced when the 


matrix order is raised. We proceed to prove that, for a given equa- 
tion, there is a lowest order to the matrices of a complete set of 
rational conjugates. 
There exists a complete set of n rational conjugate matrix solutions, each 
of order n!, of the equation 

f, (x) = x"—a,x"-!+a, 2"-?+...+(—1)"a, = 0 (3) 
over the field « of characteristic 0. 

* p. 25 (2). 


Ss 
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We assume throughout that f,,(7) has no repeated zeros. Denote 
the complex zeros by a4, a,..., x,, and the corresponding matrix 


zeros by A,, ; ee . One root of (3) is 


n 


Next 
f(x) (a A,) an —p, x" 2 Pox" 3_! al ( -a yr 1», 1 0. 
where p, = @,—a,_, A, +4,» Aj+...+(—1)*Aj, is an equation with 


2 


commutative matrix coefficients each of order n and has one rational 


matrix solution of order n(n—1), viz.: 
(5) 


Evidently A, is commutative with A, regarded as being of order 
n(n—1) by repeating* A, (n—1) times in the leading diagonal. 
Hence 
F,(%)/(w—A,)(@— Ay) 

qn-2 qua" Stiga L...+(—1)"-q,, ° 0 (6) 
has commutative’ coefficients each of order n(n—1) which are 
symmetric polynomials in A, and A,. Proceeding in this way we 
reach a complete set of n rational matrices each of order n! which 

(i) satisfy the same reduced equation f,,(2) = 0 and are therefore 
all derogatory, 

(ii) are mutually commutative. 

Moreover, in the repeated factorization of f,,(7) we may replace the 
matrices at each step by the original complex zeros ay, a,..., x, and, 
since multiplication is commutative, 

(iii) the elementary symmetric functions of the matrix roots 
A,,...,A, are the same as the corresponding elementary symmetric 
functions of the complex roots a4, a,..., x,- Hence A,, Ay,..., A, 
form a complete set of rational conjugate matricest each of order 

We shall often use the same symbol to denote matrices of different orders. 


+ i.e. they satisfy conditions (i) and (ii). 
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n‘ which are solutions of (3). Even if f,,(x) is reducible or has repeated 
zeros, the result holds provided that the coefficients are commutative. 
If the coefficients are matrices of order m, then our solutions are 
matrices of order n!m. 
Certain equations have a complete set of rational conjugate solu- 
tions in matrices of order less than n!. Thus 
(7) 
has not only 














0 —1 0 —l1 —!1 
6 #*—J i = 








obtained by the preceding method, but also the complete set of 
conjugates 

= 1] —1 —1 —l1 —] 0 
A, 100], A,=|-—1 O 1], A;s= 0 —1 —1]. (9) 

010 it 0-2 —I-—l |! 
In both (8) and (9) the elementary symmetric functions of the roots 
are also those of (7) but there is an important distinction; in (9) 
A, = Aj—2, A, 1 — A}—A,, but these relations no longer hold for 
the solutions (8). These relations are replaced by 

(A,— Aj+2)(A,+-Aj+4,—1) = 0, 


of which each factor is a divisor of zero. Next, 


The necessary and sufficient conditions that an equation of degree n 
shall be normal* is that it shall have a complete set of non-derogatory 


conjugate matrix solutions each of order n. 


* i.e. a Galois equation. 
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If f(z) = 0 is a normal equation, every complex root is a poly- 
nomial in any one, i.e. 

Tae) (vx—a){a—r,(a)}...ja—r,,(a)} 0 (10) 
where 7;(x) is a rational polynomial in the complex root «. We may 
identify this form for f,,(7) with the form (3) by means of rational 
operations and the fact that f,,(2) = 0. Hence, since there is always at 


east one matrix root A, of order n, we may replace « by A,i 1e argu- 
least t t A, of order n, yrepla by A, in theargu 


ment and f(x) -_ a—r,(A,)}...{e—r,,(A,)} = 0, (11) 


so that A,, 7,(A,),..., 7,(4,) are matrix roots which form a complete 
conjugate system. 

Conversely, if all the matrix roots are of degree n and if A, is non- 
derogatory,* any matrix commutative with A, is a polynomial in A,. 
Since the conjugates must be commutative with A, it follows that 
they are polynomials in A,, and f,,(7) may be written in the form 
(11), which, by a repetition of the first part of the argument, may 
be identified with (10). Hence f,,(x) 0 is normal. We see that 
Sylvester’s theorem that, if A is non-derogatory,+ all matrices com- 
mutative with A are polynomials in A and the theorem that in a 
normal equation all the roots are polynomials in any one of them are 
merely different ways of stating the same fact. 

If a polynomial relation 4(a, 0 connects the complex roots 
of f,,(x) 0, then some $(A,, Ag,...,A,,) is either zero or a divisor of 
zero in any matrix representation of the roots. 


For & x,), defined as the product of all the poly- 


nomials obtained from ¢ by operations of the symmetric group in the 


suffixes, is zero. Hence, since A,,..., A, is a complete set of conjugates, 


TI 4(A,.....4,,) = 0. 


Therefore some factor $(A,.,..., A.) is either zero or a divisor of zero. 


n 


If in any matrix representation 4(A,,...,A,,) = 0, then 


p. 94. The matrix A, of order n of § 2 is non-derogatory and can be trans- 
formed into the A, above. 
+ A may be non-derogatory even though f,(«) is reducible. P. Franklin, 
Ex. 3, shows that «° = 1 has the complete sets of roots 
= 0 0] ‘O —] 
( 0 —1)3 (1 o> ( l 0): 


the members of the first set are derogatory whilst those of the second set are not. 
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For [] ¢(A,.....4,) = 0 is symmetric in the roots and is 
s 


therefore unchanged when we replace A,,...,A, by a,...,«,. Thus 


the ’s being some permutation of «,,...,«,. As it is rational, we must 
also have eas 
” . ( Nyseeey Hy Ts 0. 

If A be a matrix of order n satisfying an irreducible normal equation 
f(x) = 0 of degree n, then there is no polynomial in A which is a 
divisor of zero. 

For the latent roots of any polynomial g(A) are g(a,) 

Hence by a transformation, possibly irrational, 


g(x) 


If this is zero, then every g(a,) = 0 (s = 1...n), and, since f,,(2) is 
irreducible, f,,()/g,(v) and therefore g(A) is zero and not merely a 
divisor of zero. It follows that 

The realm «(A) obtained by adjoining the matrix A of order n, solution 
of a normal equation of degree n, to the field x is a field which is simply 
isomorphic with the field (x). 

This is an extension of a theorem of Pierce* in which he proves a 
corresponding theorem for cyclic equations. 

If the representation of the roots of a normal equation of degree n is of 
order a multiple of n, then the realm «(Ay,...,A,,) may contain divisors 
of zero. 

For A, can be transformed into a direct sum of matrices each of 
which is of order and satisfies f,(~) = 0. Hence, if A, = B,(A,), 
A,, will be similarly expressed as a direct sum of matrices of order x. 
Such representations are therefore mere duplications, in the leading 
diagonal, of the normal representation and will not be regarded as 
of true order exceeding n. If the order is a true multiple of x, then 
TJ] {4,—8,(A,)} is symmetric and is rational and therefore zero. 
Hence A,—f,(A,) may be a divisor of zero. 

3. The resolvent matrix 

The general number in the realm «(A,, Ag,...,A,,), obtained by 

adjoining the conjugate roots A,, Ag,...,A,, of § 2 to the field of the 


Hn e 


* p. 25. 
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coefficients, «, is 


where x; is an indeterminate. There are n!—1 associates of V; which 
are formed from it by permuting the A’s by means of every possible 
permutation of the suffixes. It will be proved in § 4 that matrices 
T.. exist such that 

Tj; A,T,; = A;, T;;'A;T,;; = Ai 

T;;'A,T,; =A, (p #14,)). 

Hence the associates of J, being obtained from it by transformation 
satisfy the same characteristic and reduced equations which are 
therefore of degree n!, and V, is non-derogatory in the case now con- 
sidered. Further, the A’s are commutative and so are the V’s. Hence 

The characteristic equation for the resolvent V, is normal, and V, and 
its associates form a complete set of rational conjugate matrices. 

[t follows that the characteristic equation for V, is also a complete 
Galois resolvent for f,,(7) = 0. A further consequence is that J; 
and A; are all rational polynomials in Jj. There are two cases to 
consider: 

(i) if the characteristic equation of J; is irreducible 

n(A,, Ao,...,.A,) K(V3) 
is a field, and it is impossible to transform all the A’s simultaneously 
into direct sums of matrices of order less than n!; 

(ii) if the characteristic equation is reducible, a rational matrix 
T’ can be found such that 7'—'V; 7 is the direct sum of non-derogatory 
matrices of order less than n!, the characteristic equations of which 
are the irreducible factors of the characteristic equation of J. 

Also, since every J; is a rational polynomial in V,, each Vj, and 
therefore every A, is simultaneously transformed by 7’ into a direct 


sum of matrices whose orders are the same as the degrees of the 


irreducible factors of the characteristic polynomial of V,, that is of 


the complete Galois resolvent. Further, since every transform of 


a zero is also a zero and since 


then 


Hence f(A;,,) = 0 (p ‘). Also, since the A’s are commutative, 


) 


so are the A,’s, and the set A, (7 n) forms a complete set of 
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matrix solutions of f(a) == 0. As this is supposed irreducible and 
as the transform of a conjugate set is also a conjugate set, the set 
A, (i = 1,...,n) may be taken to be the same matrices arranged in 
a different order. Evidently the complete resolvent splits into r 
irreducible conjugate factors each of degree s = n!/r, but since the 
characteristic equation for V;,,, is different from that for V;,, and each 
is irreducible, it is impossible to transform J;,,—>Vj,,. Hence the 
effect of transforming by 7’ is to separate the n! conjugate matrices 
V; into classes in which the leading elements of each class form a set 
of conjugate solutions of a Galois resolvent. In Ex. 8 
aif = Ss 
—1-2 0 0 90 

] 2 1 1 
—1 0 0 —l1 —4 —% 

0 —1 0 —2 -% 

0 0 —1 0 —2 —3 








The complete resolvent is 
(a3 + 6a?+ 5a-+ 1)(a3+- 62?+5x—13) = 0, 
which is reducible although V, is non-derogatory. The other roots are 
V, = 2A,+3A,+4s;, Vz = 3A,+A,+24s, (13) 


V,=A,+34,424, = 24,+4,4+345, Ve = 34,+24,+As. 
(14) 


T, T- are respectively 


" 4-9 —5 —1 —10 —4] 
me.’ 


ae Se 


] —6 —1 
+ : y 4 3 
2 2 
] 











] 2 
4 —: 1 6] 





These transform JV, and its conjugates into T-%7, T-V,T, 
T-1,7, T-V, 7, T-V, 7, TV, 7, which are respectively equal to 


] 0 0 —!1 
wut sf 


2 1 —5 


4 0 ] 
1 -—3 —2 
—3-—3 1 





A 
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The leading elements in the transforms of Vj, V3, V3 all satisfy the re- 
solvent #°+62?+5x2+1= 0 whilst the lower elements satisfy the 
other factor of the complete resolvent, i.e. 23+ 627+ 52—13 = 0. The 
effect of transforming A,, Ay, A; by 7' is 
T-1A,T = A,+4,, T-"A.T 
where on the right-hand side we refer to the matrices (8) in their orde1 
along the leading diagonal from left to right.* 

It follows from this that there is a least value to the order of a 
complete conjugate rational matrix representation of the roots of 


an irreducible equation and that this order is a divisor of n! For this 


representation K(V,) is a field, whereas for truet representations by 


matrices of higher order, e.g. n!, «(Vj) is not a field but contains 


divisors of zero as we have noted. This fact is sometimes useful. 
In (8) and in the field K(A,) 


with resolvent 
{fV—(A+1)2H7 


Since Vis non-derogatory it follows that «(V,) in matrices of order 


* P. Franklin's Ex. 3 (p. 25) is explained by the fact that, when the charac- 
teristic equation is reducible, solutions may be constructed for which V is 
expressed as a direct sum of elements some of which satisfy the same reduced 
equation and not, as with our V,, the conjugate resolvents. Such solutions 
will be derogatory. 

t i.e. not mere repetitions of representations of lower order. 
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six contains divisors of zero and the order of the representation is 
unnecessarily high. 


4. Transforming matrices and the Galois group 

The transform of any matrix root is also a root and the transform 
of a complete set of conjugates is also a complete set of conjugates. 
Also* any two roots of the same reduced irreducible equation can be 
transformed into each other. We proceed to show that for the 
representation of order n! of the conjugate sets there exist rational 
matrices which interchange any pair of roots whilst leaving the others 
unaltered. The proof is written out for the case n = 5 but is quite 
general. 

There is a matrix of order (n—1), having as elements matrices of 
order n, which interchanges A, and Ag (of § 2). 


Write 


where p, 


and let 


where 6, is a matrix of order n = 5. 
We require that PA, = A, P and PA, = A,P. 
Now PA, = A, P requires that 
8, At = pd, A®— p23, A?+ p36, A—pydy 
which, using (15), may also be written 
A4§, = A%8, py—A*8, po + Ad, py—Sy Py. (17) 


Evidently (16) is also satisfied by 5, A, 5, A®, 5, A*. Further 55, 43, 5, 
must also satisfy (16). Again PA, = A, P requires that 
5, p,+8, = Ad), —8, po+8, = Ad,, 8, P3t+5, = Ads, 
—8, py = Ady. (18) 


* Wedderburn, Lectures on Matrices, p. 114. 
+ The solution of (16) involves only rational operations. 
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For these to be consistent 
—5, py = Ad, = —Ad, p,+A75; = —AS, p,+A?5, p.+ Ad, 
— Ad, p3+A*8, pp—A*8, p,+A%9,. 


Hence, if 5, be chosen to satisfy (16) and 6,, 53, 5, be calculated 


from (18) in terms of 5,, the P will interchange A, and Ag. 


Since the equations found by multiplying on the right by powers 
A are left linearly dependent, the equation 
pxyt—p,xA+p,xA*—p,xA?+2xA! : (19) 
is singular.* 
A particular solution is easily found if we use the fact that 


xy 


where 


and it is easily verified that sucl 
0 OO 
0 
| 
0 
0 
where s, is the sum of the nth powers of the complex roots. Hence 
P can be calculated by purely rational operations. 
There is a transforming matrix of order n! which interchanges A, 
and A, and leaves Ag,...,A, unaltered, the A’s being regarded as 


i 


matrices of order n! 4 
43 
0 


Xo 


2] 
2 





Y2 Y¥34 
where each of the elements a, 8, y interchanges A, and Ag. 


* A. R. Richardson, Proc. London. Math. Soc. (2), 28 (1929), 416. 
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Now SA, = A;S requires that 
4 G+ % = 1% —F2Bi +93 9%. Fa 43 = G1 %2—G2 Bo +s V2, 
By qi +B2 = %4, —Biq2+Bs = %, 
Mnt+V72 = Bi —7192—73 = Be, 
—%4 93 = 91 %—G2 B3 +43 Ys 
—Bi Is = % 
—7193 = Bs. 
The solution is a, = 8B,=—y3, 1% ==), —~ 8 =” == ©. 
P 
Hence S= | P | is the type of solution. 
P 
A repetition of the argument proves the result. Since «(A,) is a 
field, we may apply an induction to f,,(a)/(e—A,) and then rewrite 
the elements of the transforming matrices of order (n—1)! as scalar 
matrices of order n and so prove that matrices 7;; exist, of order n!, 


which interchange any pair of roots A;, A; whilst leaving the others 


unaltered. It is evident that the transforming matrix operators 
T;;\(_ )T;; when applied to elements of the realm «(A,, Ag,...,A,) 
form a group which is simply isomorphic with the symmetric group. 
Further, it is easy to show that any 7;; is the product of a fixed S;; by 
a number in the realm «(V,). The unit element of the group consists 
of all numbers in «(J}). 


When 7' exists which transforms every J; into a direct sum of 


matrix elements of orders in the leading diagonal, Jj,...,); may be 
chosen so as to have as leading elements a set of conjugate solutions 
of the same partial resolvent. The matrices 7;; 
J. themselves are resolved by 7' into direct sums of matrices 


which interchange 


of order s,* and the partial matrices which transform the leading 
elements of Jj,...,V, evidently form a group of transforms simply 
isomorphic with the Galois group of the equation. The relation of 
this group to the symmetric group is indicated by the way in 
which 7’ expresses the 7;;’s as direct sums of elements not all in 
the leading diagonal. 

Such matrices may themselves be chosen to form a group,} and 


* A. Voss, p. 93, not all, however, in the leading diagonal. 

+ H. Hasse, T'rans. American Math. Soc. (34), 189 and also Max Deuring, 
Algebren, 60-7, and B. L. van der Waerden, Gruppen von linearen Transforma- 
tionen, 71. 
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this leads to the interesting conclusion that the matrix roots can 
be put in a canonical form. For a group of matrices can be put into 
such a form by transformation and the corresponding matrix roots 
will then be in canonical form. 

Thus for the quadratic x2?—a, «+a, = 0 a canonical form for the 
O 


Galois group is | 
0 1 


| 0) ; ‘ 
; ) and the corresponding canonical roots 
QO — 


are 


1 1 


l(p2 ° \ 1 2 
$a, 3(A; ~~" »(@{— 4a) 
2 7 
2 yy 244 


By using (20) and P we can find a matrix which will transform 7;; 
into the corresponding canonical matrix and therefore the set of 


conjugate roots into a canonical set. We have proved that 


If the Galois group of an irreducible equation is of order s, then a 
complete set of conjugate rational matrix roots A,,..., A,,, each of order s 
can be found, by purely rational operations, which are transformed into 
each other by a group of matrices, each of order s, simply isomorphic 
with the Galois group. Further, the realm (Ay, Ag,...,A,) 18 simply 
isomorphic with the Galois field and, if the matrix representation is of 
order greater than that of the Galois group, and if V, is non-derogatory, 


then «(V,) may contain divisors of zero. 


In Ex. 7 the group of transforming matrices may be taken to 
he Land Qio, Gis, Gog: G23 Vigo respectively as 
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-—l —1 
] 

—] 0 
] : ; 0 

— | y 0 
] 0 








For the complete set of solutions (9), 7' transforms the above thus: if 


—l1-1 3 0 12 
Qies = 1 0 —2], Qiss = |—1 —1 1), (22) 
0 0 1 0 olf 


which together with the unit matrix form a group, then 


a 0 ey 
T 10,5 T | 0 ig | 7-29... bas vg | 


Qise 0 123 
Kor 4 
0 2 Qi 


1 0 


t 


T Qos T = ) TP Qj23 T 
4 Al ‘yy ! Case 0 
1 32 
pei TO ies 


If we transform (22) into a canonical form, e.g. 


010 {0 0 1)\ 
Ps 00 } Piss 10 0) 


1 0 0) ‘0 10 


then a set of solutions in canonical form is 
1) ale 
4 —6 21, 


—] 5 17 


] 


A, r 





in which the rational discriminant, 7, is put in evidence. 


5. Remarks 

The characteristic equation of the matrix g(A,,...,A,) has the com- 
plete set of conjugate roots obtained by operating on the suffixes of 
the A’s by the substitutions of the group of the matrix representation. 
This offers a simple way of expressing the equation of squared differ- 
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ences, the resolvent, the general number in the field, the discriminant 
of the field and similar numbers in determinant form. Further, an 
ideal can be expressed as a matrix and the theory of classes of ideals 
and of reciprocity may be regarded as the theory of a set of equations 
each of which has solutions which are transforming matrices for some 
of the others. Thus 2?—m = 0-and x7—n = 0 are such equations if 
m is the norm of a number in «(vn) and if 2 is the norm of a number 
in «(vm), whilst x?—m = 0 is self-involutory if the quadratic char- 


acter is-+ 1. The equation x«?— 2 = 0 is also self-involutory, for 
Vit 


the solutions given in § 2 are transformed into each other by Q)os3, 


dv 


¢132> 1. 














both of which are solutions of 2?—2 0. 

Finally, our methods will apply to unilateral equations with matrix 
coefficients even if non-commutative and also to equations having 
repeated factors. It is interesting to note too that a solution of 
vax+tbatact+d 0 is 

ba) ba 


a-} 





PHOTONS, ENERGY, AND RED-SHIFTS IN THE 
SPECTRA OF NEBULAE 
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1. Introduction 
In a recent paper Milne has constructed a complete dynamics of a 
particle, free or constrained, in the presence of the substratum or 
‘smoothed-out’ universe.* It represented a substantially new attack 
on the problem of the origin of the laws of dynamics and gravitation: 
whereas the usual procedure is to extrapolate laws of nature, known 
from observation to hold locally, so as to build up a model universe 
and so a cosmology, the new procedure is to extrapolate phenomena, 
from which alone laws of nature can be derived, fitting these pheno- 
mena together so as to construct a complete model universe and then 
determining the laws of nature by analysis of the totality of pheno- 
mena possible in this model. 

The extrapolation of phenomena to form a model universe was 
originally carried out by means of the ‘cosmological principle’.+ But 
it has recently been shown that the elementary properties of systems 


satisfying the cosmological principle can also be derived by use of 


the ‘sample principle’, according to which it is only necessary to 
ensure that local phenomena are everywhere the same; world-wide 
pictures are then automatically identical.t It is indeed immaterial 
whether the actual universe obeys the cosmological principle; all we 
want to know in the long run are the laws of nature holding locally, 
and we can attack this question by deriving them as the condition 
that local phenomena are in effect capable of being fitted together 
to form a seJ]f-consistent universe. If the laws of nature so derived 
did not agree with local experience, then we should have a proof that 
the whole actual universe did not obey the cosmological principle. 
The dynamics constructed by this method does, however, agree with 

* Milne, Proc. Roy. Soc. 154A (1936), 22-52. I wish to acknowledge the 
considerable assistance which I have received from Professor Milne in the 
preparation of this note for publication. 

+ Milne, Relativity, Gravitation, and World-Structure (Oxford, 1935), Chap. 3: 
referred to here as WS. 

t Whitrow, Zeits. fiir Astrophys. 12 (1936), 47-55. See also WS § 134; 
Hubble, Halley Lecture (1934), 4. 
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the dynamics observed to hold good locally in the sense that it 
reduces to the Newton-Einstein dynamics near the observer. But in 
the new dynamics energy is an invariant, not the time-component of 
a 4-vector, and so is the same for all fundamental observers. 

The purpose of the present note is to extend this dynamics to 
photons, by the usual process of regarding a photon as the limit of 
a particle of zero rest-mass moving with the velocity of light. The 
result is to give, without any appeal to the quantum theory of 
radiation, the fundamental relation* E = hy. But the astronomical 
setting of the result is novel. It exhibits the red-shift of the nebulae 


alternatively either as a Doppler effect or as a gravitational red-shift 


due to loss of energy as a photon, in moving from a nebula to 
the observer, runs down the apparent potential-gradient of the sub- 
stratum. Either interpretation is equally valid, so that the ‘redden- 
ing due to age’, to use the common phrase, can be considered at will 
as kinematical or dynamical. It will be understood that in speaking 
here of a gravitational red-shift we are referring to the red-shift due 
to the field of the substratum or general distribution of nebulae; the 
general theory of relativity (not here appealed to) would predict a 
possible additional red-shift, due to. the fact that we are near one 
condensation (our own galactic centre) and that the photon has 
originated near another, so that there may be an extra (positive or 
negative) shift due to the consequent additional difference of gravi- 
tational potential. 


2. The energy of a particle 

Let O be an observer associated with a fundamental particle, that 
is with a particle of the substratum. If O adopts Euclidean space, of 
coordinates x, y, z with himself as origin, and Newtonian time f¢, 
measured from the natural zero of time, then in his description, the 
substratum has a distribution of particle-density ndxzdydz, where n, 
the number of particles per unit volume near (x,y,z) at epoch f¢, 
according to O, is given by 

a Bt 
c3(¢2— P?/c?) 

* Kermack, McCrea, and Whittaker, Proc. Roy. Soc. Edinburgh, 53 (1932), 
31-47, and Synge, Quart. J. of Math. (Oxford), 6 (1935), 199-204 have proved 
that when a photon is emitted by one atom and absorbed by another in a 
gravitational field the ratio energy/frequency is the same for emission and 
absorption. 


(1) 
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P being the vector (x, y, z), B a constant of zero physical dimensions, 
and c the speed of light. The velocity V of any particle P of the 
system, at epoch ¢ according to O, is given by 


V= Pi. (2) 


In the presence of this substratum, any free particle moving 
through position P at epoch ¢ with velocity V in O’s description 
undergoes an acceleration dV/dt given by 


dV Y 
dt ~~ (P—Vd)> G(é), (3) 


where X= #—P4J/c2, Y= 1—V%/c?, Z = t—P-V/c? (4) 


‘Se (t—P-V/c?)? (5) 
~ XY (#—P?/c2)(1—V2/c2) ~ 


In his new dynamics Milne takes G(é) = —1 and assigns to any 
particle of rest-mass m moving with an arbitrary velocity V through 
the position P at epoch ¢ (all relative to O) the energy E = mc*£! and 
inertial mass M = m&!. For P = 0, i.e. when O is chosen to be at the 
fundamental particle close to the given particle, ZH reduces to 
mc?(1—V?/c?)-+, M to m(1—V?/c?)-*, which are Einstein’s expres- 
sions. In all cases EH = Mc? as in Einstein’s dynamics. 

It is readily seen, however, that, if the system admits the concept 
of energy as an invariant for all fundamental observers, then, whatever 
the value of G(€), the energy of a particle must be given by Milne’s 
formula. Kor the only two independent functions of P, V, ¢ which 
are invariants are €, X and hence the energy of a particle must be 
of the form ¢(€,X). Consequently the energy of a slow-moving 
particle in the neighbourhood of O is given by ¢(€, X) where £ ~ 1, 
X ~ #. It follows that, if our definition of energy is to be compatible 
with the well-established principle of conservation of energy for 
slow-moving particles in the neighbourhood of our own galactic 
centre, then ¢(€,X) must be independent of X. Finally, for any 
particle passing O the energy must be given by Einstein’s formula* 
mc?(1—V?/c?)-#. Since the value of € for P = 0 is (1—V?/c?)-1 we 
find that ¢(€, X) must be of the form mc?é?. 





and é 


* Einstein deduced this formula dynamically, whereas the present con- 
siderations are purely kinematical. However, if the definition of energy here 
adopted is to be of physical significance, it is clear that for a particle passing 
O the expression must be of the same algebraic form as Einstein’s. 

3695.7 Y 
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3. The energy of a photon 

We consider a photon as a particle moving with the velocity of 
light. Milne has shown, from purely kinematical considerations, 
that such a particle must necessarily contain in its description a 
parameter, of the dimensions of a frequency, additional to the 
specification of its position and direction of motion.* Our object 
is to link this parameter with the energy of the particle, by purely 
kinematic arguments, i.e. without assuming the quantum theory 
of radiation. 

We first observe that for a particle moving with the speed c the 
invariant € is infinite. Hence, since the equation of motion (3) can 
be written in the form 

dV Z G(é) 
a (P—Vt)—, 
dt a” £ 
a photon may be regarded as a free particle describing a rectilinear 
trajectory with the velocity of light if and only if 


GE) _ 9, (7) 





, (6) 


lim 

+00 

Consider now a photon passing through the fundamental observer 

O. Let v be O’s measure of its frequency. The photon is to be 

regarded as the limit of a particle whose speed |V| +c and whose 

rest-mass m—> 0. In order that # may have a finite limit, the ratio 

m/(1—V?/c?)! must have a finite limit. This limit must depend on the 
description of the photon,t i.e. on v. We put accordingly 


m 

= lim : 8 

uv) nm (1—V2/e2)! (8) 
|\Vi-—ec 


The energy EH of the photon, being an invariant, must have the 
same value for the given photon at its given position and epoch for 
all fundamental observers. A particular fundamental observer is that 
one past whom the photon is flying, and we can calculate E by deter- 
mining its value as calculated by O. Accordingly 

2 ln2 
, mc?(t—P-V /c? 
E = lim ( Ic) 
m—>0 ( P=0 
=e 
.2 
= c*u(v). (9) 

* Milne, Zeits. fiir Astrophys. 6 (1933), 83. 

+ In effect we are assuming that the energy of a photon is completely 
specified by its frequency. 





#—P?/c2)i(1—V/c2)} 
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Now let a photon be emitted at time ¢, reckoned by O, from a 
nebula or fundamental observer O’ in the direction of O. Let its 
energy at O’ be E’; let its frequency to O’ be v’; and let its frequency 
to O be Vv. Then by (9) E’ — c2ulv’). (10) 


But the same value Z’ must be assigned to the energy by O, since 
the energy is aninvariant. Putting OO’ = P, we see that O calculates 


the energy as 2(¢_P-Vie2 
bal me?(t—P,-V/c?) (11) 
kad (#?— P§/c?)*(1—V?/c?)* 
Vie 


In this expression, since the direction of V is that of the vector —P,, 
P. 


V= —c-__®. 12 
Bi (7) 





Combining (11) and (12) we have 


2 i IP,|/c 
EF = aad li mec 0 
~ no (1—V2/e2)! (@—P2/e2)h 


lVi-e 


aes, t+-|Py!/c ‘ 


But P, = Vjt, where V, is the velocity of O’. Hence 


ee 





Equating (10) and (14) we have 


ur’) = ply ieee i. (15) 


But by the ordinary kinematic theory of the Doppler effect, estab- 
lished by light-signalling analysis only without recourse to dynamics, 
the frequency v observed by O is less than that, v’, observed by the 


(16) 


receding O’ in the ratio 1—|V,|/c 
i "| 1 + ¥, I Eel 


By (15) and (16) ov) (17) 


, 
v 


Since O is arbitrary the ratio must be the same for all fundamental 
particles. Call it h/c?, where h is a universal constant.* Then 

pv) = hv/ c, 
and so by (9) E=bh. (18) 


* By assuming that for a photon m/(1—V?/c*) — pu(v) we have necessarily 
found that / is a secular invariant. The possibility of h varying with the epoch 
has been discussed by Chalmers and Chalmers, Phil. Mag. 19 (1934), 436-46. 
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We are thus led immediately to the isolation of the quantum con- 


stant h. 


4. Red-shifts in the spectra of nebulae 

The energy of the photon when at O’ is EZ’ = hy’ and when at O 

its energy is H = hv. Thus 

EK’ uae v’ =a 1+ |V,j/e 

z->7 lime 
Consequently the energy of the photon diminishes as it passes from 
O’ to O. Its frequency, either to O’ or to O, remains constant during 
this passage, but the invariant energy decreases from E’ to E. It 
advances with constant speed and so we may say that its potential 
energy in the gravitational field of the substratum diminishes as it 
passes from O’ to O, for Milne has shown that in his revised dynamics 
the invariant expression F counts in both what we are disposed to 
call kinetic energy and what we are disposed to call gravitational 
potential energy in the field of the substratum. Thus we have at will 
either a kinematical or dynamical description of the red-shift. 

It must be noticed that in the present formulation the quantum 
relation H = hv is not an invariant relation. ZH is an invariant, the 
same for all observers at any given event in the photon’s history; v is 
different (in general) for every different observer. E is constant for all 
observers for a given photon at an assigned position; v is constant for 
a given photon for an assigned observer. The equality ZH = hv is just 
satisfied for each observer as the photon passes that observer. 

It may be repeated that the history of a photon as here analysed 
is that of a particle moving with the speed of light in the smoothed- 
out model. We have thus taken into account only the general 
gravitational field of the universe. The effect of condensations would 


to 


(19) 


require a separate treatment. 
5. Summary 

The revised dynamics recently derived by Milne from purely 
kinematical considerations is applied to photons. The relation 
E = hv is shown to follow, by considering the passage of a photon 
from a nebula to the observer in the gravitational field of the 
smoothed-out universe. It is further shown that the red-shift may be 
considered either as a kinematic Doppler effect, or as a loss of energy 
consequent on passage through the gravitational difference of poten- 
tial between nebula and observer. 
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1. In the theory of canonical matrices two matrix pencils 
I, = pA,+oB, and I, = pA,+cB,, 


each with m rows and » columns, are said to be equivalent if two 
constant non-singular matrices P and Q of degrees m and n exist 
och ies 6 PLO=f, UP #101 #@. (1) 

The fundamental result is that two pencils can be reduced to the 
same canonical form, if and only if they are equivalent. Suppose 
now that [, and J, are equivalent; then it is possible to solve (1) for 
P and Q, and the natural question arises: What is the most general 
solution of this equation? It is easily seen that this problem is 
equivalent to finding the most general pair of matrices P, Q which 
transforms a given pencil [= pA-+-oB into itself, i.e., which satisfies 
the ti 

1e equation Pre =f, (2) 
or, comparing coefficients of p and o, 

PAQ= 4; (2a) 
PBQ = B. (2b) 

In this case we shall say that the pair of matrices (P,Q) is an 
automorphic transformation of I’, and we propose to determine all 
such automorphic transformations; in particular, we shall express 
the number of linearly independent ones among them in terms of 
the invariants of I’, i.e., in terms of the invariant factors and the 
Kronecker indices of I’. 

We may make some convenient assumptions regarding 

[' = pA+oB, 
without restricting the generality of the investigations. 

First, we shall assume that A is of the same rank as TI itself, a 
condition which can always be fulfilled by a linear transformation 
of the variables p, o. For, let po, 79 be such that 

A, = pA+oB 
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is of maximum rank, i.e., of the same rank as [ itself, and put 
B, = p, A+o, B, 
where p, and a, are only subject to the condition 
Po Pil 49g 
Op 2% 
If we now introduce new variables p’, o’ by the transformation 
, , 
P = pop’ +pie 
yy , 
o = 09p +90, 
we see that I can be written as 
, , 
I’ = p’A,-+0’B,, 
and that it now has the property required. Again, the automorphic 
transformations remain the same, since the equations 
2 — 
PA,Q = A), (4a) 
2 “ 
PB,Q = B, (4b) 
are equivalent to (2a) and (2b) on account of (3). After this pre- 
liminary remark we shall replace the two homogeneous variables 


p, o by one variable, A, and write the pencil in the form 


lr = AA+ B, 
where A is of the same rank as I’. 
Next, may be replaced by any pencil I, which is equivalent to [’. 
For, let L=srr (|S| 40; |7| £0). 
If (P, Q) is an automorphic transformation of I’, then 
(SPS-1, T“1QT) 
is an automorphic transformation of [,, and vice versa. Thus a 
(1,1) correspondence is established between the automorphic trans- 
formations of [ and [,. In particular, we may assume that I, is the 
canonical form* of I, which we write in the form 


a. 2 
Ny = 1,+2,+4;, 
Pp 
where = ps (ajAm,) (0 << m, < mg <... << m,), 


t 
Tr, = A,+- M,, 
qd 
nM ." a - a 
Yr — psy (B;A,,,) (0 ~ Ny = No 2s coe SN Mq)s 
Jj 
using the ‘direct sum’ of matrices which is defined as follows: 


* See the footnote on the next page. 
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4+ht+...4+0 =z I, = diag(P,, J..... 5) = r, : 


and, when [| = [, =... = [. = I, we write 


r4+r+...40 = (1). 
In (5.1) and (5.3), A,,, is the typical singular sub-matrix* correspond- 
ing to a row vector of minimal degree m; annihilating I’, e.g. 
- 
A, = K! Ag={f1l Alj,..., 
1 


and in general d 
1 


A 
1 (6) 


s+1,s 
We remark that A, has s+-1 rows and s columns. The pencil I, in 
(5.2) is the ‘non-singular core of I’ of type kx k, say, and there is no 
loss of generality in assuming that the matrix coefiicient in [, is 
the unit matrix; for since that coefficient is of the same rank as J}, it 
must be non-singular and may be removed as a matrix factor. 
The Kronecker minimal indices are exhibited in the canonical 
form, viz: 
a, times m,, x, times mzy,... for the first set, 
B, times n,, B, times ng,... for the second set. 
These numbers together with the invariant factors of [, are the in- 
variants of I: they are the same for all pencils equivalent to ’. In 
what follows we shall assume that I is already in canonical form, 
ie. LU = f 
Since !Q| 4 0, we may put R = Q-! and write (2) as 
PY=PR (|P| #0; |R| #0) (7) 
or by (5) P(L.+h+h) = (.+1h+h)z. (8) 
We now partition P and R in conformity with the three isolated 
sub-matrices of I’, i.e. we put 
P = [P,), R= [R,;] (4,f = 1, 2,3). 


* See Turnbull and Aitken, Canonical Matrices (Glasgow 1932), 121; and 
W. Ledermann, Proc. Edinburgh Math. Soc. (2) 4 (1934), 92-105. 
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Hence (8) resolves into nine partial equations 
| 
Z PT; = 0, Ry, 
or, in full, 
ae = ode om. re = 
iA L=Tky» iv R= Re, (vii) Ast, = 1 Ry, 
@PF5=-R2, ~e&GE=—h2 we) F1,=—Tee! © 
(ili) Py T= Ty Rg, (vi) Peel, = Ty Rye, (IX) Pog Ty = Ty Regs, 
and it is the object of the following pages to give a complete solution 
to these equations, the m?+-n? elements of P and R being regarded as 
the unknowns. Only equation (v) seems to have received attention 
in the literature. For substituting for [, from (5.2) and suppressing 
unnecessary indices we can write this equation as 
P(AI+M) = (AI+-M)R, 
whence r= & 
and Pa — MR. 
Hence Pe =— MP. 
The solution of (v), therefore, involves the finding of the most general 
matrix P that commutes with a given matrix M. 

This problem was first solved by Frobenius and has since been 
treated by several authors.* It was found that the number of 
linearly independent solutions of (v) is 

t, = t = e,+3e,+5e,+7e,+..., (10) 
where e, is the degree in A of the vth invariant factor of the matrix 
AI,,+M,. In what follows we shall obtain similar results for the 
remaining eight equations of (9) to which we shall refer later simply 
by Roman numerals. The total number ¢ of parameters in the general 
solution of (7) is equal to the sum of the sub-totals ¢,, ¢g,..., t, giving 
the numbers of parameters in the solutions of those nine equations 
(9). It will be found that three of these numbers are zero, i.e., that 
the corresponding equations have only the trivial solution in which 
all unknowns vanish; and the final result will be 


9 
TT Py =r p a; w,(m;—m,;+ 1) + 2, PiBilns—m +1) + 
+k(> a+ > Bi)+ Dd x; B(m,+-m,). 
t J tJ 


* Frobenius, Berliner Sitzungsb. (1910), 3-15, where other references may 
be found; also D. E. Rutherford, Proc. Ak. Wet. Amsterdam, 35 (1932), 870-5. 
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2. Before solving the commutantal equations § 1 (9), we shall make 
some remarks on the typical singular sub-matrix (§ 1 (6)): 


1 A 
1 


A, = (1) 


where J, is the unit matrix of degree s, and the dot below or above it 
indicates a row of zeros. 
First, we wish to determine the row vectors and column vectors 
that annihilate A,. Let 
= t= {€,, &5,---» €,} 


be a column vector with s elements and 


Y = [No M+ Ns] 
a row vector with s+1 elements. We then prove the following 
lemma: 

Lemma 1. The equation 

A,z=0 (2) 
admits only of the solution x = 0, i.e. the columns of A, are linearly 
independent; the general solution of 

yA, = 0 (3) 
is y = d(A)[1, —A, (—A)?,..., (—A)*], 
where (A) is an arbitrary function of and 

u, = [1, —A, (—A)?...., (—A)*] 
is a vector of degree s. 

Proof. The rank of A, is s, because, on cancelling the first row of 
A,, we obtain a minor of order s identically equal to unity. In the 
set of homogeneous equations (2), the number of unknowns therefore 
equals the rank, and by the fundamental theorem of linear equations 
the system has only the trivial solution x = 0. In the set of equations 
(3), however, the number of unknowns exceeds the rank by unity 
and the most general solution is a scalar multiple of any particular 
solution which may be taken to be the vector u, defined above, since 

u,A, = 0, (4) 


as is easily verified. 
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Lemma 2. If P is a constant matrix with s+-1 rows, then the equation 


a P= 0 (5) 
is impossible unless P=0 


The proof follows immediately on writing out (5) in full and equating 
the coefficients of 1, A, A?,..., A® to zero. 

3. We now turn to the discussion of equation § 1 (9). Consider 
the simple case PA, = A, R, (1) 
which is solved by the following theorem. 

THEOREM 1. When s’ < 8, equation (1) has only the trivial solution 


a, R= 06; 2) 


when s' > 8 and s'—s = d (say), then the general solution is 


P= [bij les 1,s+1 (: 


and R = [$:-jle.s 

where do, $4,---, bq are d+1 arbitrary constants and 
ima 3-9 

pa \(i—j < 0). 

For instance in the case s’ = 5, s = 3, d = 2, we have 
i‘ |, Youee 
be Me Re 3 ae 
$s 1 po ¥ ai : dy 

$2 $i bo 
n . 2 1 
en a 7 Se 
L- + Ss ; 
Proof. On premultiplying (1) by u, we obtain by § 2 (4) 
(u, P)A, = 0, 








i.e. the vector u, P annihilates A,. Hence by Lemma 1 (§ 2) we have 


8 

u, P = g(A)[1, —A, (—A)?,..., (—A)?]. (5) 

The elements of the vector on the left-hand side are polynomials of 

degree not higher than s’. Hence, comparing the first elements of 
either side, we see that 4(A) is also such a polynomial. 
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Now, when s’ is less than s, equation (5) is obviously impossible 
unless (A) is equal to zero and hence 
u, P = 0, (6) 
which, by Lemma 2 (§ 2), implies 
P = 0. 
Equation (1) then becomes 0 = A, R which entails 
R= 0, 
since the columns of A, are linearly independent. This proves the 
first part of the theorem. 
Next, when s’ > s, (5) can evidently be solved and ¢(A) will, in 
general, be a polynomial of degree d = s’—s; 
$(A) = $o—biAt+h2r*+...+-ba(—A)% 
say. Putting 
P=([p,] (¢ = 9,1.,...,8’; j = 0, 1....,8) 
and comparing the jth elements of either side of (5), we obtain 
Poi t+Pry(—A)+Paj(—A)?+... ADei(—A)” = bo(—Al +-64(—AP + 
+¢4,.(—A)i#?+...4+46(—A)*4 (7 = 0,1, 2,..., 8). 


Hence 
Poj —— Prij =— we Pj-1,3 = 0, 


Pi3 = Fo Pps = $1, +) DPjraj = Par 
Pjsa+1,j = + = Pes = 9; 
0 (t—j < 90), 
Pis = (Fi (0S 1-j < A), 
0 (—j> 4), 


which proves the statement (3) regarding P. In order to determine 
R we substitute § 2 (1) in (1): 


I, &) 
Pial’*|+| > |) =fal*]+]- |lz. 
AT La =e 
Comparing the constant terms we get 
Pi j=], 
id= [7 )8 
which, on premultiplication by [. Z,], becomes 


[. £)P[;] = 2. 


8 
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The matrix R is thus expressed in terms of P, and it is readily seen 
that R is obtained by cancelling the first row and column of P. This 
proves equation (4). It is easily verified that (4) and (5) are also 
sufficient conditions that P and R should satisfy (1). 


Corollary. When s’' = s, the general solution of 
fa, = A,B 
is P= $9 4,41; R= $f, 
involving one parameter dy. The solution is non-singular, if and only 
if do #~ 9. 
We now come to the equation 


P(aA,) ~ (a’A,)R, (7) 


where (aA,) = A,+A,+...+A, (a times repeated). Writing P and 
Ras ‘ 
DP ees = 
P = [Py] [a (' 
) = 
we see that (7) resolves into aa’ matrix equations 
PyA, = Ay By 
of the type which we have just considered. Hence, if s’ < s, (8), and 
therefore (7), is impossible, and, if s’ >s, each equation (8) has 
s’—s-+1 linearly independent solutions, and the number of para- 
meters in the general solution (7) is consequently 
aa’ (s’—s+1). 
In particular, when s’ = s and a’ = a, P,; and R;; must be of the 
form 


Py=¢ubn By=¢ul, (69 = 1,20), 


ij 


where the ¢,; are «* constants (Corollary to Theorem 1). Introducing 
a matrix 
® = [44], 


we can write the result in this case as 


P = [$4 Jor] = OX Jay 

R = [$i] =OxI, 
using the familiar notation for Zehfuss (or Kronecker) matrices. The 
solutions P and @ are non-singular, if and only if ® is non-singular. 
This enables us at last to obtain the complete solution of (i) 


(§ 1 (9)): 
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THEOREM 2. T'he general solution of 


PLE; (i And] = [3; And] 2 
is of the form 
Py Ry 


P on Py Pr iP us Ry Ry. 


ee ... ‘. = 


pl 
involving in all 
g t, = fe, eh — My +2) (10) 
: 
arbitrary parameters. 


The proof is obvious from the preceding results. For, again, we 
partition P and R in a suitable way so as to resolve (9) into the 


equations P; (cx; Ani) on hee Ay.) Rij, 


which have already been discussed. When i < j, we have § 1 (5.1) 
and therefore P;; = 0; R,; = 0. Further, P,; and R;; are of the form 
Py, = ©; X Dy i433 Ry =O; T,; and P and R are non-singular if and 
only if ,, ,,...,0,, are non-singular. 

The solution of (ix) (§ 1 (9)), namely 


qd. ; Ns qd ; 
PLS: BsAn)] = [2 BAn)]R, 
does not involve any new difficulties, since the transposition of this 


equation leads us back to the previous case. The number of para- 
meters in the general solution of (ix) is therefore 


tb = > Bi Bnj—n;, +1). (11) 
joi 
We shall now consider (iv), viz. PI, = [, R, 


ie. Pl, = [: (a;Ap,)|R (12) 


where T, = A, +m, 
is a non-singular pencil of type kxk. As before, (12) resolves into 
a number of equations of the kind 
PT, = A,, BR; (i = 1,2,...,5 a) (13) 
j 
where s; stands for m, (x, times repeated), m, («2 times repeated), etc. 
It is easy to see that (13) is impossible, because premultiplying it by 


U,, we get (u, PT, = 0. 


sit 
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Since [, is a non-singular pencil, it follows that u,, P; = 0, and hence, 
by Lemma 2 (§ 2), P; = 0. Equation (13) now becomes 
Oo=A, R,, 
which entails R,; = 0, 
the columns of A,, being independent. The solution of (iv) therefore 
contributes no parameters, i.e. 
0. (14) 
We get a different result, however, when [, and [, are interchanged, 
as is the case in (ii), viz. 
PT, = [ B, 
i.e. Pi x (;Am,)| = (AI,+M,)R, 
i 
which may be resolved into > a; equations of this kind: 
> > | L, . \ | F 
P,A,, = PAV | + = (AU,,+M,)R; (i 
\ > r% J . 
where s; takes the values m,, mg,..., m,, the m; being repeated «; 
times. Hence comparing coefficients of A we get 


[|= | 
P| | = M,R, 


where, for simplicity, we have dropped the suffixes of the matrices 
P and R. Let 


(16) 


R = [1y, Fe)... %,] 
introducing column vectors pp, 7,.-., 71, %2;--- for the columns of P 


and R each having k rows. It is easy to see that postmultiplication by 


ye or F has the effect of cancelling the last or the first row. 


81 
Hence equation (16) becomes 
[Por Prs-+s Psa) = [Pa T25+++5 %ei)> 
[P1; Pos++s Ps) = [M11 Mj, 12,---; My 19,1; 
and by eliminating 1, 79,..., Tg, 


[2 1» Pos-++s Pe] = [Mh Po, Mj, P1;---, My Po,-1]) 
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whence 
Pr = MyPo P2 = MP, = MEM;---» Py = MiDo,--- 
(v = 0, 1...., 8; 
T= Po, T2 = M,Po,---> Ty = MEM, (uv = 1,2...., 8%), 
where the k elements of the vector p, remain arbitrary. We have 
therefore: 

THEOREM 3. The general solution of the equation 

PA, = (AL.+M,)R; 
P; = [Po My Po, ME Po,--- Mi Po], 
R; = [Po My Po ME Do:---, Mi'Po) 
involving k arbitrary parameters, namely the elements of pg. 

Since (15), or (ii), is equivalent to } «; equations of type (15’), we 
have at once = kD xy. (17) 
As before, we see that (vi) (§ 1 (9)), ie 
and (viii), ie. PI,;=TR 
are merely different forms of (ii) and (iv) to which they are reduced 
by transposition. [, is then replaced by a pencil of the same type as 
I (§ 1 (5.1) and (5.3)) and the order of the factors is reversed, while 
I, and 1; play the same role since the only property we have used 
was that [, was a non-singular pencil of type kx k. 

By analogy, we obtain therefore: 


qd 
t= kD Bs (18) 
t, = 0. (19) 
We shall now show that (vii), i.e. 
PT, = £& 2B, 
ie. PLS: BAs,)] = [E- @%Amd]B (20) 
J v 
is satisfied only in the trivial case P = 0, R = 0. For (20) splits up 


into a number of equations 
Pi An, = Am, Bij (21) 


mi 


is of the form 


which after oe ae by u,, yield 
Um; Py = —_ 0, 


whence Um, Pi = 
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the rows of Aj,, being linearly independent (§ 2, Lemma 2). Hence, 
as before, it follows that 
P;=0 and R,; = 0. 
We therefore have t = 0. 
It remains to solve (iii), which we write as 


P[S: (%:Am)| = [E- BAn)|2 


v J 


and which reduces to ( > a;)( > f;) partial equations 
\ 4 7\ j 


PsA 
each of them occurring a;8; times. Substituting for A,,, and Aj,, from 


nj 
§ 2 (1), we get 


Pl ™|+[7 |} = Old Re 


mi 


ao AL Buy, (24) 


m 


= [Z, “| Rj; 


“j 


Puv - 


Pu,v YS Tutt 
From these two equations we infer 
Puv+i —_ Putty 
or, replacing v by v—1 and iterating the equation, 
Puy = Pu+i,v-1 ben eee Puth,v—h — veo Pu+v—2 Say, 
depends only on the sum of the suffixes, and 


ie. the value of p,, 


Vv 


similarly 


Tuy — Tutl,v-1 a ge: 0 Puv = Pu+y-2 
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We have therefore proved the following theorem: 
TueoreM 4. The general solution of P,;A,,, = Aj, Ri; ts of the form 


T do py : ’ , Pn, 7 


P. . oy dy ° . = dm; +1 


y 


es Pn; ° fo 
} bo $1 . e E Pini—1 7 
1 2 . . . Pm, 


R., 


Wy 








Ldn; Pnj +t 4 : . Pangy~4. nj+1,m; 
involving m,+-n; parameters $q, $45---» Pm,+nj-1- 

The number of parameters that occur in the general solution of 
(iii) is therefore 


= ye ;Bj(m,-+-n,;) (: . " ei (29) 


a > Speen, 


Ry now add up a sub-totals t,, ty,..., fy as given in § 1 (10) and 
(10), (11), (14), (17), (18), (19), (22), and (29), and we find that 
Hh number of linearly independent solutions of 
: PY =TR (30) 
is equal to 
€ = 64-4... Fo, 
= t+k>a;+k ZB Z maj —my +1) + 


+ 2 BiBilrs—mit 1)+ p> a; B,(m;+n,;). 


Moreover, a method has been given for actually obtaining all 
matrices P and R satisfying (30). Since 
t, = t, = ts = 0, 
it follows that in the scheme § 1 (9) 
P,=0, Rp = 0; Ps =0,R3= 0; Pyy=0; Ry = 0, 
so that the general solution is of the form 
BP, te ‘ Ry ‘ 
P= |Pu Pr | R= |Ry Ry .« |. 
Py P52 Py Rs Rp Rs 
Further, P and R are non-singular if and only if the matrices 
P.,, Poo, Pyg and R,,, Ry, R3g are non-singular, and we have already 
found the necessary and sufficient conditions that those matrices 


should be non-singular (see p. 285). 
3695.7 U 








PATH-GEOMETRY AND COSMOGONY 
By D. D. KOSAMBI (Poona) 
[Received 2 June 1936; in revised form 31 August 1936] 
THE equations proposed by E. A. Milne (1,2) for the movement of 


free particles in his expanding universe may be written in the form 


da*\? 
. 
d2x7 4 F dx? > (za) : 
iene ("aah SHO = © (1) 


= wt dx’ r 

(: >: iia) 
: — /dxi\2 
ior Benh— D(a] 


These are obtained by putting 2° = ct in Milne’s original equations. 





(i = 1,2, 8). 


The equations given above may be taken to give the paths of a three- 
dimensional space and discussed as such by methods outlined else- 
where (3,4,5,6), but I shall replace them by an equivalent four- 
dimensional system in which x° appears as an additional coordinate. 
We change to a new parameter 7 by means of the identities 
d d ld 1 d? d? i # Pd 
cdt dx® dr c2dt? d(x)? (@°)2 dr? (#°)8 dr 
In these formulae, as also in all those which succeed, the dots over a 
letter indicate differentiation with respect to the parameter 7. The 
equations (1) now become 
gt yo - P pout — > (xt /°)2 G(é) 0 (i 
<_< £4 aii . F ~~ 9 7 — 2 l 
(a°)? (#°)8 4g (x°)?— > (x*)? 
These are multiplied by (#°)?, and may then be written as 
) X 
where Y = (#°)?— > (#*)?, X = (x°)?— > (x*)?, 
9) and X being equivalent respectively to c?Y and c?X in Milne’s 
notation. To simplify further, let the tensor g;; have the components 


x x° 


g_@ ét—(xt— =a Vee) =0 (i =1,2,3), (3) 


Joo = 1, Gui = -i, Ju = Ji0 = 0, Fi. 0 (7 # §). (4) 


Take p‘ to be a four-component contravariant vector which, in the 
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particular system of coordinates chosen, has the value x‘. Then, 


Pi = GirP’> Ly = Yin’, 3= pt, =p, = Gis %'p', (5) 
inj oe asag 2 
X=pp,=gyr'r, Y=HE=gytt, §= is 


The indices now and henceforth range over the values 0, 1, 2, 3, and 
the tensor-summation convention is to be followed throughout for 
repeated indices; but the use of the tensor g;; for the purpose of 
raising and lowering indices is to be confined to the two vectors 
p’ and # only. 

We have only three equations in (3); the fourth, for #°, remains to 
be specified so as to give the simplest possible system. This is 
obviously effected by choosing 

24 
a(?—e@| =0, =p (6) 

Setting apart the case # = 0, which corresponds to an instanta- 
neous change in the three-dimensional universe, we can now combine 
all four equations as follows: 

Zitat = 0, 
(i = 0,1, 2, 3). (I) 


Sane 


X 


Using the notation and formulae given in my former papers, the 
differential invariants of the path-space defined by (I) are given 
for the kinematic case G(€) = —1 by 


PY Ge) 


a 


. Oat , 
ites _ OC ik = 9, 
P| = —XY3i+ Ypip, + ¥i#,— 3p", 

X? Ri, = X(8j,.4;—3) ,,) +-p'(p; %,—p, &;), 

XH Rip) = XB}. Ij —8} Jar) +P" (D5 Ir — Pr Gn) (7) 
They may, of course, be calculated for all G. The question of the 
existence of a metric is of particular interest and has, to a certain 
extent, been settled by Walker. But his metric W(£,X)VvY cannot 
be valid even if the integral for his function @ (8,484) be given a 
meaning by some sort of a limiting process. I proceed to show 
that for no metric can the equations (I) represent geodesics when 
G(é) = —1. 
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In fact, if such a metric exists, it must be homogeneous in # (6) 
and equivalent to a Finsler metric or to one homogeneous of degree 
zero. But such a metric f is to be obtained as a solution of the system 
of partial differential equations 


of, +#f,=0,  fe—dolehr = 0. (8) 
These have the conditions of integrability 
Pif:, — 0, Rijfer == @.... (9) 


From the homogeneity of f, we have 2’f,, = kf. Substituting this in 
the identity p’P;f,; = 0, it follows that 
if = 3 Pf, (10) 
But the restriction p/ Ri, f., = 0 gives us 
P'S. Ps = ¥f.5- (11) 

These may be combined to give f.; = kfp;/3, which has the unique 
solution f = A(3)*. This, however, is a trivial solution, since the 
condition for regularity |f,,.;| A 0 is not fulfilled, i.e., the equations 
(1) cannot represent the geodesics for such a metric. 

The general case G(é) 4 —1 may also be discussed by similar 
methods, though this will be of no particular value. It may be 
pointed out, however, that f being any homogeneous metric re- 
ducible to one of Finsler’s, or of null degree, any arbitrary ¢(f) is also 
a possible metric. In addition, there exists along the paths (or, 
for the Finsler space, must be assumed) the integral f = constant. 
If any metric homogeneous of degree zero exists, it is not included in 
those deduced by Walker. The space G(é) = —1 is purely affine, 
having no metric. 

The interest of the equations (I) lies primarily in the fact that they 
lend themselves to generalization, and are tensor-invariant. As a 
matter of fact, equations derived from any ‘least action’ principle 
have the advantage of being tensor-invariant, this being from a 
formal point of view the most important use of the calculus of 
variations. Combination of (I) with the De Donder-Einstein-Mayer 
equations may be brought about by adding the terms Ij, #/%*, the 
“, being Christoffel symbols associated with a general tensor g;; 
which is to replace that in (4). For the Maxwell tensor, the terms 
— F‘¢* must also be added to a‘. The vector p' may be replaced by 
one of a more general type. In all this there must be the proviso 
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that a system of coordinates exists—if not simultaneously for the 
whole space, at least at every point—such that the tensor g;; in them 
may be expressed in the form (4), the vector p‘ having at the same 
time the components 2’. 

It is not clear, however, whether the original field equations given 
by De Donder, Einstein, and Mayer for g;;, Fj, and the curvature 
tensor formed from the I, should be modified, or whether their 
equivalent may be derived from equations in terms of the new 
tensors calculated for the modified a‘. Even when a metric does not 
exist, it might still be possible to construct a purely affine field 
theory (7). In all this, it must be kept in mind that there are two 
distinct spaces associated with the equations (I): the space of which 
(I) represents the paths, and the space of which g;; is the funda- 
mental tensor. 
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/ ON THE LEBESGUE CLASS OF INTEGRAL 
FUNCTIONS ALONG STRAIGHT LINES 
ISSUING FROM THE ORIGIN 


By V. GANAPATHY IYER (Madras) 
[Received 25 May 1936] 


1. Introduction. In an interesting paper* on the mean values of 
functions regular in a strip, Hardy, Ingham, and Pélya have proved 


that, if the integral 
«o 


F(a) = [ \fle-riy)| dy 
—w 

converges for x = a and x = B while f(x+-ty), regular in the strip 
x < «x <8, satisfies the condition 

f(x+iy) = Ofexp(e*')} (0 <te< = ) (1) 
uniformly in the strip, then J(x) converges for each x in a < « < B 
and J(x) < max[J(«), J(8)]; it is also shown that log J(x) is a convex 
function of x. 

1.1. In this paper similar results are proved for an angley under a 
less restrictive condition (C) (§ 1.3) than (1). In §2 the results are 
proved for restricted values of p by using the properties of analytic 
functions alone. In §3 these results are extended to any p > 0 by 
using the theory of sub-harmonic functions. The methods of this 
paper can be applied to prove the convexity property as well, though 
this point is not discussed in the body of the paper. 


* Hardy, Ingham, and Pélya, Proc. Royal Soc. A, 113(1926,7),554, Theorem 7. 
+ Ihave not been able to prove, by the methods of this paper, the analogous 


results for a strip without an additional assumption that f | f(w+-iy)|? dy is 


-— & 
uniformly bounded in a < # < B for some fixed 7’. This does not involve 
any loss of generality under condition (1) as is shown in the paper quoted in 
the footnote above. But I have not been able to show that this is the case under 
the less restrictive condition (C) (§ 1.3). 
t The following line of argument enables us to prove the convexity pro- 
x 


perty. Let H(0) = f | f(te"®)|” dt. Consider the function H,(z) of § 3.2, Lemma 2. 


0 
When € varies in (0, ©), {H,(z)* represents, under the conditions of Theorem 3, 
a uniformly-bounded normal family of sub-harmoniec functions. The function 
v(z) contemplated in § 3.1 (iv) is precisely H(@) as defined above. Hence H(@) 
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1.2. Definition and notation. In what follows, we mean, by a line, 
a half-line issuing from the origin. We mean, by an angle, one or 
other of the two infinite regions enclosed by two half-lines issuing 
from the origin; when we wish to speak of the magnitude of the angle 
we consider that angle which maps out the region under considera- 
tion. Let f(z) be a function regular in the interior and on the sides 
of an angle A. We shall denote by M(r,A,f) the maximum of |f(z)| in 
the sector cut out by A on the circular region |z| < r. 


1.3. We require the following Lemma:* 
Lemma 1. Let f(z) be a function regular in an angle A of magnitude 
a. Let* lo. M A 
og MAD < o, © 


a 


lim 
Let \f(z)| < A on the sides of the angle. Then |f(z)| < A throughout 
the angle. 

2. THEorEM 1. Let f(z) bea function regular in an angle A of magni- 
tude x. Let it satisfy the condition (C) in A. Let f(z) belong to the class 
L,, (p @ positive integer) along the sides of A. Then it belongs to L, 
uniformly along every line in A. 

Proof. Let 25 = r,e be any fixed point in A. Define P(r) for 
0<r< rT, by the relation 

P(r) = [sgn f(re’)|-” 
where, as usual, 
f(2)| [sen f2)] = fle) (fe) 4 9) 
and sgn f(z) = 0 (f(z) = 0). 
Let z = re‘ and 


1 r 
Mid = [ [f(eA)}” P(r A) dd = ef | [sep P("e1 dt. 


“ 


Now H(z) is a function regular in A. By hypothesis it is bounded 


is sub-harmonic and constant on the curves @ = constant. Therefore by (v) 
§ 3.1, H(8) is a convex function of 6. By applying this result to z~*!? f(z) (a > 0), 
we conclude the possibility that e®H(@) be convex in 6 when « > 0. Hence 
by another result proved by Montel in the paper quoted later we conclude that 
log H(@) is a convex function of 0. 

* R. and F. Nevanlinna, Ac. Soc. Fennicae, 50. 

+ I am indebted to a referee for replacing my original slightly less general 
and cumbersome conditions by the elegant hypothesis (C) in the Lemma and 
the succeeding theorems. 
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along the sides of A. Moreover, if f(z) satisfies condition (C), so does 
H(z). Hence, by Lemma (1), H(z) is uniformly bounded in A. But, 


at 2 = 2p, 


|H(z9)| = | \f(te)| at. 
0 


Since z, is any point in A, the theorem follows. 

2.1. THeoreM 2. Let T be a system consisting of a finite number of 
lines. Let « be the magnitude of the greatest of the angles between con- 
secutive members of the system. Let f(z) be an integral function which 
satisfies (C) in each angle A between consecutive members of 1. Let 
f(z) belong to some class L,, (p > 1) on the lines of 1’, where p might be 
different on the different lines. Then f(z) = 0. 

Proof. Let « be any given positive number. Consider the function 

¢ re 
=z | f(zA) dA = ec” | fi(te') dt, 
0 0 
where z = re’’. G.(z) is an integral function and satisfies condition 
(C), when f(z) does so, in each angle A between consecutive members 
of I’. Using Schwarz’s inequality, we easily find that, whatever 
p > 1 be on the sides of A, |@,(z)/z| is bounded along these sides 
as |z|—>0o. Hence by an easy modification of Lemma 1, we con- 
clude that |@,(z)/z| is bounded in each angle A between consecutive 
members of I. Therefore |G,(z)/z| is bounded in the whole plane as 
z|—0. Therefore G(z) is a polynomial of degree one at most. 
Differentiating with respect to « and putting « = 1, we find that 
zf(z) is of the form az+h. Since f(z) is an integral function it must 
reduce to a constant which is zero since f(z) belongs to some L,, along 
lines of [. 

2.2. The following corollaries are easily derived from Theorem 2. 

Corotuary 1. If the function f(z) be of order p less than m/x or if it 
be of order w/x and minimal type, it will be identically zero under the 
hypothesis of Theorem 2. 

CoroLuaRy 2. If an integral function of order one and minimal type 
belongs to L,, (p > 1), along a WHOLE line, then it is identically zero.* 

For, we can take the origin at some point on the line so that we 
get two half-lines inclined at an angle z. So we can apply Corollary 1. 

* When p = 2, the result of Corollary 2 has been proved by Paley and 
Wiener, American Math. Soc. Collog. Pub. xix, 13, Theorem XI. 
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CoroLuaRy 3. If an integral function of order one belongs to some 
L,, (p > 1) on two INTERSECTING WHOLE LINES in the plane, then it is 
identically zero. 

Here, we can take the point of intersection of the two lines as the 
origin. The greatest angle between the consecutive lines is less than 
a. Therefore we can apply Corollary 1. 


3. We shall use the theory of sub-harmonic functions to extend the 
above results to the case p > 0. We require the following from the 
theory of sub-harmonic functions:* 


3.1. A real continuous function u(z) of z defined in a domain D is 
said to be ‘sub-harmonic in D’ when the value of the function at any 
point z, in the interior of D does not exceed the mean value of the 
function along all circles of sufficiently small radius round 2, ; 

(i) a sub-harmonic function does not attain its upper bound in 
the interior of D unless it reduces to a constant; 

(ii) if w and v are sub-harmonic in D, so are u-+v and cu, where c is 
a positive constant; 

(iii) if {w,,(z)} is a uniformly convergent sequence of sub-harmonic 
functions in D, then the limit function u(z) is sub-harmonic in D; 

(iv) let {uw} be a uniformly-bounded normal family of sub-har- 
monic functions, and v(z) be the function equal to the greatest limit 
of the set {u(z)} at z where wu varies over the functions of the family, 
then v(z) is sub-harmonic in D; 

(v) let w(z) be a harmonic function, and v(z) be a sub-harmonic 
function taking constant values on the curves u(z) = constant, then 
v, considered as a function of uw, is a convex function of wu; 

(vi) let f(z) be regular in D, then |f(z)|? (p > 0) is sub-harmonic 
in D. 


3.2. We next prove the following lemmas: 


Lemma 2. Let ¢« be a fixed positive number. Let f(z) be a function 
regular in an angle A. Let z = re‘ and 
H(z) =r | |fed)|” dd = | |f(te)|r dt (p > 0). 
0 


€ 
0 


Then H(z) is sub-harmonic in A. Let g(z) be regular in A. Then 
(eH (z)| is sub-harmonic in A. 


* For these properties, see P. Montel, J. de Math. (New Series), 7 (1928). 
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Proof. We have 
n 
\e!@H_(z)| = lim = 


no =< 


Dp 


zlpeg@lp (= | 
n 








uniformly in any finite part of A. Hence, by (iii) and (vi) of § 3.1, 
|e%H_(z)| is sub-harmonic. The first part of the lemma follows on 
taking g(z) = 0. 

3.3. Lemma 3. Let f(z) be regular in an angle A of magnitude « and 
satisfy condition (C). Let H(z) be defined as in Lemma 2. If 
H(z) < A on the sides of A, then H(z) < A in the whole of A. 

Proof. By Lemma 2, |e?H-(z)| is sub-harmonic; therefore* 
U = log H(z) is sub-harmonic. Let B be any contour lying in A. 
Let V(z) be the function harmonic in the interior of B and equal to 
U(z) on B. Let G be the Green’s function belonging to B. Then 


V(z) = 2. [ Vis ) eas 


27 | on 
B 
1 , ¢ iG 
=~ | U(s) —ds, 
2a . on 
B 


where ds is an element of the arc of B. Now U(z)—V(z) is sub- 
harmonic in the interior of B by (ii) of § 3.1, since —V(z) is sub- 
harmonic. Also U(z)—V(z) = 0 on B. Hence U(z)—V(z) < 0 by 
(i) of § 3.1 in the interior of B. Hence for z in the interior of B 
1 E G 
U(2)< [ U(s) 2" as. 
2a | on 
B 
From this point, the proof now follows that of Lemma 1. 
THEOREM 3. Theorem 1 is true for any positive p. 
Proof. We consider the auxiliary functions defined in Lemma 2. 
By hypothesis the conditions of Lemma 3 are satisfied by H,(z), and 
H(z) < A, some constant, on the sides of A. Hence putting « = 1, 


€ 


the theorem follows from Lemma 3. 

4.1. THrorEM 4. Theorem 2 remains true when p > 0, provided it 
is the same p for all the lines of T. 

Proof. As before, we consider the auxiliary function defined in 
Lemma 2. By Lemma 3 applied to each of the regions between 


* See P. Montel, loc. cit. 40 (footnote 1). 
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consecutive lines of I’, we conclude that H(z) is bounded in the whole 
plane. We, now, require the following 

4.2. Lemma 4. A sub-harmonic function bounded in the whole plane 
reduces to a constant. 

Proof. Let u(z) be the sub-harmonic function. Consider the family 
{u(ze")} of sub-harmonic functions, where 0 < @ < 2x. In any circle 
'z| <r, the family is bounded and normal. Hence, by (iv) of §3.1, the 
function equal to the greatest limit of the set {u(ze)} (0 < @ < 27) 
at the point z is sub-harmonic. But this is evidently M(r), where 
z = re” and M(r) denotes the maximum of u(z) on |z| = r. By (i) of 
§ 3.1, M(r) is a non-decreasing function of r. By (v) of § 3.1, M(r) is 
a convex function of logr. By using the convexity-property of 
M(r) it is easy to show that M(r) is a constant since it is bounded. 
Hence u(0) = M(r). We can now apply (i) of § 3.1 and conclude that 
u(z) is a constant. 

4.3. Resuming the proof of Theorem 4, we find that H(z) is a 
constant by Lemma 4. Differentiating with respect to « and letting 
« = 1, we get {z||f(z)|” is a constant which can only be zero since 
f(z) is an integral function. This completes the proof of Theorem 4. 

4.4. Conclusion. It might be asked whether an integral function 
can belong to some L,, along every line. The answer is in the affirma- 
tive. We know* that there is a function g(z) tending to zero along 
every line. The function 


g(2)+9(—2)—29(0) 
22 


belongs to L,, (p > 4) along every line. But all the same, no func- 
tion can belong to some L,, uniformly along every line unless it is 
identically zero. For, the auxiliary function H(z) of Lemma 2 would 
be bounded in the whole plane so that by Lemma 3 it reduces to a 
constant. Hence, as before, f(z) = 0. 


* See P. Dienes, Taylor Series (Oxford, 1931), Chap. X, 346. 








NOTE ON A THEOREM OF HELMHOLTZ 
By J. HODGKINSON (Oxford) 
[Received 20 April; in revised form 28 September 1936] 

A FLUID occupying infinite space is at rest at infinity. Within a 
closed surface S there are an expansion 6 and vorticities (€, 7, €); 
6, €, n, € are everywhere zero outside 8S. Then the velocity com- 
ponents (w, v, w) at a point (x, y, z) are given by Helmholtz’s formulae 

oo eH a 

Ox ' Oy a2’ 


oF oH 


> 


+= 


ti o da'dy'dz’, 
+ r 


“Pg 
{| | > dx'dy'dz’, 
~~ « r 


| | i dx'dy'dz’, 
> 


= | [ | 3 du'dy'dz' ; 
et aN a 


the integrations are taken throughout the space enclosed by S, and 
r is the distance between the point (#,y,z) and the point (a’, y’, 2’) 
at which the volume-element of the integral is situate.* 

It has been suggested that, although there is no doubt that any 
fluid motion may be represented by the formulae, some restrictions 
must be placed on the values of 0, €, », in order that the formulae 
may represent a fluid motion. This note provides a discussion of 
this point. 

The formulae may be interpreted in two ways. If they are required 
to represent a fluid motion for all time, i.e. 0, £, , ¢ are functions 
of the time ¢ as well as of the coordinates, it seems that the formulae 
are nothing but an alternative presentation of the fundamental 


* H. von Helmholtz, J. fiir Math. 55 (1858), 25-55. The notation used is 
that of H. Lamb, Hydrodynamics (6th ed. 1932), Chap. VII, 208-10. 
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differential equations of hydrodynamics. But we may regard the 
formulae as representing a fluid motion for one instant only, say 
the initial moment ¢ = 0. With this interpretation 6, &, n, ¢ are 
functions of the coordinates only. For precision of statement we 
replace the surface S by two surfaces S,, S,. Within S, there is no 
volume throughout which @ vanishes; within S, there is no volume 
throughout which €, », ¢ all vanish. I shall suppose that 0, and €, », f 
are regular functions within S, and S,, and, of course, that 
og as 


n+gte = an 


I shall also suppose that the fluid pressure p is a function of the 
density p, and that the body forces may be derived from a potential. 
I shall show that the only further restriction necessary in order that 
Helmholtz’s formulae may represent a fluid motion at an ‘ordinary 
moment’ is that the value of @ must tend to zero at all points on 
S,, when these points are approached from within S,. [By the phrase 
‘ordinary moment’ I mean a moment at which the acceleration at 
a point of the fluid is everywhere continuous. | 

A very brief discussion of the first interpretation of the formulae 
suffices. Within S,, S, the values of 0, £, », ¢ are those of a fluid 
motion. They must satisfy the relations* 


1 Dp 
p Dt 
mA.) fou , nou, Cou 
Dt\p pcx pey péz 


5 +0 = 0 (A) 


Di\p p ox p oy p oz 
D(t\ _ € ew , » ow ¢ ow 
Dt\p} — -p dx p oy | pez 


m2) = E ov , 1 Wv ¢ ov 


But, since 

(€,,¢) = 4 curl(u, v, w), 6 = div(u, v, w), 
these may be regarded as three differential relations satisfied by 
p, u, v, w in any fluid motion. [The equations (B) are not indepen- 
dent. If they are differentiated with respect to x, y, z respectively, 
the equation found by adding the differentiated relations is satisfied 


* E. J. Nansen, Messenger of Math. 3 (1874), 120; Lamb, loc. cit. 205. 
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identically in virtue of (A) and the relation ~ 


0€ on , % ex, 


Ox Oy @& 








They have the same kind of dependence as the equations éz/éx = P, 
éz/éy = Q, where 0Q/éx = @P/éy.| It seems therefore that in order 
to find 6, €, n, ¢ so that Helmholtz’s formulae may represent a fluid 
motion, we must determine a fluid motion to which @, &, n, ¢ belong. 
In other words, we have merely restated the analytical problem pre- 
sented by the differential equations of fluid motion. 

We turn to the second interpretation of the formulae. The in- 
tegrals ©, F, G, H are of the same form as the integral defining the 
gravitational potential of a volume-distribution of mass, and the first 
derivatives of such a potential function are continuous everywhere. 





Hence wu, v, w are continuous everywhere. Thus the motion defined 
by the formulae is kinematically possible. This is the sense in which 
I interpret some sentences of Lamb’s account.* 

We cannot, however, ignore dynamical considerations altogether. 
The differential equations of fluid motion are formed with the 
assumption that the pressure-gradient is continuous,t and conse- 
quently so are the accelerations of the particles of fluid. It is con- 
ceivable that, at a ‘singular moment’ of the motion, the acceleration 
might become discontinuous at the points of some surface within the 
fluid without the collapse of the motion. It might well be true that 
Helmholtz’s formulae would represent a motion at a singular moment 
unless the values of 6, , n, € were restricted in some way. 

It is clear that, if we show the accelerations of the particles of fluid 
to be continuous at points on S,, S,, we need consider them at no 
other points. To do this we shall require to make use of the theorem 
concerning the discontinuity of the normal component of the attrac- 
tion of a surface-distribution of matter, and the corresponding 
theorem concerning the discontinuity of the second derivatives of 
the potential of a volume-distribution.{ [I state the latter, which 
is less familiar.] Suppose the density of the distribution within a 
surface S is p. Let V, be the potential at points outside S, and V; 

* Loc. cit. 204, 206. 

+ In forming the equations (B) we assume the pressure-gradient differen- 
tiable. 

t G. Kirchhoff, Vorlesungen iiber mathematische Physik: Mechanik (dritte 
Auflage 1883), 178. 
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the potential at points inside S, and let 1, m, n be the direction- 
cosines of the outward-drawn normal at a point P on S. Then, at P, 


aM, ON _ 

i tie 
Me OM snip. 
Oxdy Oxdy 

It will be found necessary to assume the existence of initial values 
of 06/ét, 0&/Ot, On/et, 0f/et at points within S,, S,. The justification 
of these assumptions is postponed. 

We consider the x-component of acceleration, i.e. 


ou ou ou ou 
—+u—+v—+uw—, 
ot ¥ Ox * oy = Oz 


© oH o@ 

Ox Oy 2 

For the moment we suppose that 0, €, n, { are functions of ¢ as well 
as of the coordinates. Then 


2 — = LS | If cl . (; dx'dy'dz’ + ee If v, v=(;) d8,. 
oxot 4 ot ox\r 4n Oxu\r 


The second term of the right-hand member represents that part of 
the rate of variation of 0®/éx which is due to the motion of the 
surface S,. This surface will not move with the fluid, in general, and 
v, is written as the velocity of the surface normal to itself. Assuming 
the existence of an initial value for 00/et, we may now put t = 0. 
The first term in this expression for é°®/éxét can be represented as 
the a-component of attraction of a volume-distribution of matter. 
This component is continuous at points on S,. The second term is 
similarly represented as the 2-component of the attraction of a 
surface-distribution. The discontinuity of the normal component is 
—v,0. Hence the discontinuity of the z-component is —lv, 0, i.e. we 
may write 


where “= — 


OU, Ui _ py 
ot Oot ati 


We next turn to éu/éx, éu/éy, @u/éz. The contributions of 6 to 
these quantities are 
eet ao eo 
x2’ éxdy’ xd’ 


which we represent as second derivatives of the potential of a 
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volume-distribution. By means of the theorem I have quoted we 
find that 


Cu, OU; ___jap Ou, OU: np 


ce 


= —Iné. 


Up, OU; 

Ox 80x Cy oy dz (Oz 

Hence the total discontinuity in the x-components of the accelera- 

tions of the particles of fluid 2t points on the surface S, is 
U(v,—lu—mv— nw). 

It was noted that, in general, the surface S, does not move with 
the fluid, i.e. v,—lu—mv—nw does not vanish. In order, then, that 
the accelerations may be continuous at points on S,, it is necessary 
that the value of @ shall tend to zero as we approach such points 
from within §,. 

We consider the possible discontinuities of the accelerations at 
points on S, in the same manner. These are calculated from the 
integrals F, G, H. The important difference in this calculation is 
due to the fact that S, moves with the fluid, and the factor corre- 
sponding to v,—-lu—mv—nw above vanishes. Hence the presence 
of the vortices causes no discontinuity at points on the surface S,,. 
In this calculation we have assumed the existence of initial values 
of a€/ét, n/t, oL/at. 

We have now to justify the assumption of initial values of 0€/ét, 
én /ot, OC /at, 00/at. 

/ 

The equation ail: 

Dt\p 


Edu , nou, Feu 
W pox poy ped 
may be written in the form 


yet ate 


——U- v- 
oz ox oy 


ou é é 0 

a ~ et we 
ot Oz 
Helmholtz’s formulae give the values of wu, v, w at points within S,, 
and we are thus able to assign an initial value to é¢/ét without further 
information. 

ou ov , Ow 

Further @g= —+—+4+—, 
Ox Cy 
au av 


so that r ee a as 
wot oyot =z 


also 


ou ou 
whence —v— ) : 


De 
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To justify the assumption of an initial value of 06/ét at points within 
S,;, we must assume that the pressure p is a twice-differentiable 
function of p, and that the potential of the impressed forces and the 
initial density of the fluid are twice-differentiable functions of the 
coordinates. The same conditions are in fact required to establish 
the dynamical equation 

Dé 

Dip 
These latter conditions, however, do not apply to the prescribed 
values of 6, &, n, ¢. They operate as limitations upon the system 
whose motion is described by Helmholtz’s formulae. 

It is now established that, provided that 6 tends to zero at points 
on S,, the approach to these points being made from within §,, 
Helmholtz’s formulae represent a motion of a fluid at an ordinary 
moment in the sense that the accelerations are everywhere con- 
tinuous. If @ does not satisfy this provision, the formulae represent 
a fluid motion at a moment which is ‘singular’ in the sense that the 
accelerations at points on S, are momentarily discontinuous. 


ete 


_ € eu 
~ péx' p by’ poz 








HYPERGEOMETRIC PARTIAL DIFFERENTIAL 
EQUATIONS (II) 


By T. W. CHAUNDY (Ozford) 
[Received 1 September 1936] 


1. Euler’s equation 
In a previous paper under this title (3) I discussed the solution of the 
partial differential equation 

f(8)F(8')V = xyg(8)G(s')V (3.8 =o, v5] 

Cx oy 
and its extension to fields of more independent variables, obtaining 
their ‘general solutions’ in terms of ‘fundamental solutions’ expressed 
as multiple power-series of hypergeometric type.* In this new paper, 
presupposing the same structure for the general solution (which is 
essentially Riemann’s), I obtain an alternative expression for the 
fundamental solution which is simpler in practice and has some 
advantages from the theoretical point of view. 
I begin with Euler’s equation 


oOoT7 
py 


veal 
(c—y) —_ +n———m> (1) 
Oxoy Yy 


of which the solution is known:} or rather with the variant 
88'V = xy(d—m)(8’—n)V (2) 
obtained from (1) by writing 2”V for V and x-! for x (or, of course, by 
the similar transformation in y). The adjoint of (2), in the operators 
6,0, m8 88'U = xy(8+-m+1)(8’-+n+U. (3) 
The fundamental solution U(X, Y; x,y) is a solution of (3) such 
thatt 
(i) on 2 = X, it satisfies 8’U = Xy(d’+n+1)U, ie. it is of the 
form U = C,(1—Xy)-*; 
(ii) on y= Y, it satisfies 83U = xY(8+m-+1)U, ie. it is of the 
form U = C,(1—axY)-™*); 
(iii) ata = X, y = Y, it reduces to (1I—XY)-1. 
290 (8); or more precisely 293, § 3. 


* 3, 
+ 4, chap. 3, 54-70; chap. 4, 81-4. 
t 3, 294 (21), (23). 
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The last condition fixes the constants C,, C,, and we can state the 
conditions more compactly thus: 


U = (1—XY)"/(1—Xy)"*"* onzx= X, (4) 
U = (1—XY)"/(l—aY)"""  ony=Y. (5) 
I now suppose m, » positive integers and write 
U = (8 +1)(8'+2)...(8’ +n) W, 
so that (3) becomes 
5’(8’+-1)...(8’ +-n){b W —ay(+-m-+1)W} = 0. 
A possible solution is given by 
5W—ay(d+m+1)W = 0, 
i.e. by W = d(y)(1—ay)-—™, 
where ¢ is an arbitrary function. Thus a possible solution of (3) is 
= (8'+1)(8' +-2)...(8’-+-m){b(y)(1—axy) 9} 
= gals 1)...(8’ —n+ I){y"$(y)(1—ay)-“+?} 


y"$(y)(1—ay)-™, 


Changing the arbitrary function ¢, write this 


~ By” 


a. , = My) —zy pot 


Then, on y = Y, U = ¥(Y)(1—aY)-™», 
Thus (5) is satisfied if we take 
oy) = (1—Xy)”. 


Let us try, then, 


a on re n 
U(X,¥s xy) = (YO —xym—ayyml, 


On x = X it reduces to 
on ((y—Y)" 


sien dy" | n! 


(—Xy)| 
oie a tf 
ve mt Fra) nia—Xy))’ 

where P,_,(y) is a polynomial of degree n—1. Hence, onx = X, 
U = (1—XY)"/(1—Xy)"*+, 


which satisfies (4). Thus (6) is the required fundamental solution. 
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Working in x instead of y, we should similarly have obtained 


~ {(z@—X)" y_eyym(1—ay) m+n, (7) 
éx™\ m! J 


U(X, Y;2z,y) = 
Finally, we can show by similar arguments that 
ees om gm+n -—X ym es sé n e 
iL. Yi ege 2 (ee (8) 
dx™oy” | m!n\(1—ay) | 
is a solution of (3) and satisfies the conditions (4), (5). It is not 
difficult, of course, to verify that (6), (7), (8) are the same function. 
We can identify (6) with the known fundamental solution* 
U = (1—ay)-™+ n+1)( l —Xy)"( l —2xY)" x 


(~—X)(y—Y) 
arora (®) 





x P| —m, —n;1; 


where F is the elementary hypergeometric function, if we notice 
from (6) that (l—ay)"U is, by Taylor’s formula, the coefficient of 
h” in the expansion (in ascending powers of h) of 
(1—X[y+h(1—ay)] mfy+h(l1—ay)—Y}" 

{1—a[y+h(1—ay) ]}"*1 


a {h(w—X)+(1—Xy)(1—ha)}"{(y—Y)(1—ha) +h( _2Y)y" 
7 (l—ay)"+1(1—ha)m+t —o ‘ 





Thus (l—ay)”+"+1U is the coefficient of h” in 


>> (m!s)(n!t)(~a—X)*(y—Y )'(1—X y)"-8(1—aY )"“h" +81 —ha)-8-1, 
s=0 t=0 


where (m!s), (n!t) denote binomial coefficients. Now the terms in 
which s < ¢ are polynomials in h of degree n—1 and so do not con- 
tain h”: nor do the terms in which s > ¢. Thus, at length, 

U = (l—ay)- +") Y (m! t)(n! t)(a—X)(y—Y (1 —Xy)™-“(1—2Y)" 


t=0 
—_ (l—ay)-™" +D(1—Xy)™(1—aY)" x 
(«—X)(y—Y) 
F\ —m, —n; 1; E 
x mM, m1; Q—xyi—2¥) 
which is (9). 


* 4, 83 (22) gives the fundamental solution of Euler’s equation, which is 
equivalent, by elementary transformations, to (9). 
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2. The fundamental solution as a contour integral 
These forms for the fundamental solution become more significant 

when expressed as contour integrals. We may write (6) as 
(1—Xt)y"(t—Y)" dt 
(atm yy? 

Cc 


2miU (X,Y; x,y) = 





(10) 


where C is a simple contour including t = y but excluding ¢ = 2". 
If we write t = 7—!, this becomes 
(r—X)™(1—7tY)" dr 


2miU = 
“ (ray ry)" 





(11) 


” 


where C’ is a simple contour including 7 = x but excluding 7 = y—.* 
This is evidently the form we get by converting (7) into a contour 
integral: the symmetry of (10) in x’s and y’s is therefore clear. 

We now remark that 


u = (1l—azy)-™+), v = (l—azy)™ 
are the solutions of the adjoint pair of ordinary differential equations 
(in x) du = xy(5+m-+ 1)u, dv = xy(b—m)v; (12) 
while u = (l—2zy)-™*», v = (l—azy)" 
are similarly the solutions of the adjoint pair (in y) 

8’u = xy(d’+n+1)u, d’v = ay(d’—n)v. 
If we write these four solutions (in order) as 
E(ay), E(xy), n(xy), F(xy), 

the fundamental solution (10) can be written as 


2riV(X,¥; 2.9) = [ RoE XOm YOANN, —— (4) 
Cc 


where C is a contour including t = 0,y and excluding ¢t = 2~1,00. 
The arbitrary constants of integration which multiply the solutions 
&,,... in (14) are presumed to have been suitably chosen. 
Passing on now to the general equation in the field of two inde- 
pendent variables 
LE)F(S)V = xyg(8)G(0')\V (15) 
and its adjoint equation 
f(—8)F(—8'U = xyg(—8—1)G(—8'—1)0, (16) 


* Provided that C includes the origin and excludes ‘infinity’. 
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we construct, on the analogy of (14) and with the same contour C, 
the function 
= ; en _ dt 
2niU (X,Y; x,y) = | [> S(xE(X)}{ ¥ nly/ta(V/)}=, (17) 
r 8 
é 
where {€,(a)}, {&,(x)} are the sets of solutions (suitably paired and with 
suitable constant multipliers) of the adjoint pairs of ordinary equa- 
tions in x f(—8)u = xg(—8—1)u, f(d)v = 2 g(8)v; (18) 
and {n,(y)}, {7,(y)} are similarly chosen sets of solutions of the adjoint 
pair of equations in y 
F(—s’)ju = y G(—8’—1)u, F(8'v = y G(8’)o. (19) 
Now it is known* that, if &,,..., €,, are a linearly independent set of 
solutions of a linear ordinary differential equation 
(ap 5"+ a, "-1+-...4-a4,)v = 0, 
then a linearly independent set of solutions of the adjoint equation 
is given by é (8, 1)"-*1é,,..., €, | 
a — AF —n ! 
—_ | eg eee Gal 
the symbols denoting Wronskians formed with 5. If the solutions 
&. are chosen in this way, we have the relations 
> é,H%,= 0 (h= 0,1,...,n—2), 
E EFH4E, = og?, 
and the n reciprocal relations 
2 ESE = 0 (h = 0,1,...,n—2) \ (20) 
ry E.8"-1€, = (—)"“lap? j 
If now the &,, &,, ,, 7, of (17) are associated in the same way, the 
relation (20) enables us to show that the fundamental solution U 
has the proper values on each characteristic « = X and y = Y, and 





, etc., 


at their intersection z= X,y= Y. 
If p, g, P, Q are the orders of the operators f(3), g(3), F(8’), G(8’), 
these proper values are as follows:} 


ona = X, s’'U =0 (h<p-—1), (21) 
F(—8')?"U = 0 (p> 4), (22) 
{F(—8')—Xy G(—3’—1)}s°1U = 0 (p=4q), (23) 


with corresponding values on y = Y; 
* See, for instance, 4, 115 (12); or, in effect, 2, 382 (25). The change from the 
standard operator D to 6 presents no difficulty. t 3, 294 (19)-(23). 
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atz = Z,y = fF, 


(—8)?-(—8’)P20 = {1 (w>qor P>Q), (24) 
\—X¥)" (p=q;P=@Q). (25) 
Now on «= X we have, by differentiation across the sign of 


integration, 


2rid'U = | {> HE (até (XY n.(y/thi(¥ De 
Cc 


= 0, by (20), if kh < p—1: this gives (21); 
and, again by (20), if h = p—1, 


P 1 dt 
2n7i(—8)?-1U = | —— nAY /t)}—, 
i(—8) | xx Mula 0} 


where a (x) is the leading coefficient in (18). Thus, if p = q, so that 
aX9(X) = 1—z, 


2mil(—8)U],_~ = | LZ ally xp 
Cc 


oa [ > (nya) (7 = ¢-), 


o 
wv 


c 
i.e. [(—8)?7U],_~ = & 1(Xy)i(XY). 
Since, by definition, every 7,(Xy) is a solution of 
{F(—38’)—Xy G(—8’—1)}u = 0, (26) 
we see that (23) is satisfied. 

If p > q, the leading term a(x) is now 1, which can be obtained 
by applying the limit 2-0 in the previous value. This limit in 
(26) gives [(—5)”-1U],_. to be a solution of F(—8’)w = 0, so that 
(22) is satisfied. 

Similarly, at the intersection x = X, y = Y, if B)(y) is the leading 
coefficient in (19), 


9-74(—§)P-1(—_§’)P- =_ we 
BiB wat fa(X*)B(Y 


dt ; 

— os aa - 

| aces ifp=q,P=Q, 
é 

so that (—8)?-(—38’)?41U = (1—XY)-?: which is (25). If p > q or 

P > Q, we get the appropriate result by taking one of the limits 

x —> 0 or y > 0, which then gives (24). 
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We have still to show that U itself is a solution of the adjoint (16). 
Now, as regards the arguments x, y only, U is sufficiently written 


U = | kenny, 
: 


where u = € and u = 7 are solutions of the u-equations in (18) and 
(19). Thus, by differentiation across the sign of integration, 
f—ayR(—8U = [ (—aeao} (8 ny/0} 4 
g 
a 4) at 
= [ enyinto(—a.eenya—8yntying F 


é 
= xy g(—s)G(—3’)U. 


This gives U to be a solution of (16) and completes the proof that 
(17) is the fundamental solution of the adjoint pair of equations 
(15), (16). 
3. Comparison of forms for the fundamental solution 

I say ‘the’ fundamental solution because a theorem can be proved 
that two solutions U(X,Y; x,y) of (16) are identical if they simul- 
taneously obey all the ‘characteristic’ conditions (21)—(25). For this 
reason we can identify the expression of U as a contour integral 
given by (17) above with the expression of U as a multiple power- 
series given in (3). Let us distinguish the two forms as U, (the 
series) and U, (the integral). The formal identification of U, and 
U, can be effected by using results due to Burchnall (1). I modify 
his notation to accord with that of (3) and I use 8 not D as the 
fundamental operator. Then, as in (1), we get for typical solutions 


of (18), (19) 
E(x) = a4 s a, x*, 
s=0 


ny) =y-4 X Any®, ay) = 94 & Asy’, 


where a, A are typical zeros of (5), (5), and 


8 


§ 1 
g(a—t—1) g(a-+-t) 
- te. 
~ - Tp ~ LT iererty = 
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Thus, not forgetting the factors pu in (6), (7) of (1) we write 


2niU, = J >> (Eee toa iay td dt 





(a) F(A) t 


Applying the calculus of residues we expand the expression in { }, 
formally, in powers of ¢, retaining only the absolute term. Since a 
typical index of t is r+s— R—S, we write 

r+s = R+S =n, say, 


and sum ‘diagonally’, in r, s, and again in R, S. This gives 


5 ) ‘t i eo era 


By (27) we can rewrite 
TI g(a+s— — 
a n 


Tf fats—0 


, etc., 


where []|’ indicates that the factor given by ¢ = s is f’(a) instead of 
f(a). This completes the identification of U, and Uj. 
The form U, is free from some of the limitations imposed on U, in 


(3). In the first place U, is applicable to the ‘logarithmic’ case in 
which zeros of f or F are equal or differ by an integer. In U, the 
series must be reconstructed with logarithmic terms, but U, retains 
its form provided, of course, that the proper expressions (now in- 
volving logarithms) are put for the solutions of the associated ordinary 
equations (18), (19). 
Again, we can adapt the form U, to equations of the type 
SS)FO WV = xg(5)G(8’)V, (28) 
in which the factor y has disappeared on the right. We could deduce 
the appropriate form of U, by a limiting process, writing y” for y and 
then taking the limit h > 0: but we can equally well apply the general 
arguments of § 2 to (28) as it stands, one pair of the associated 
ordinary equations being now 
F(—8')u = « G(—8’)u, F(8’)v = « G(8’)v, (29) 
which are satisfied by powers of y with exponents involving x. 
As an example consider the adjoint pair 
58’V = 2(8+a)(5+5)V, 88’'U = a(8—a+1)(8’—b)U, (30) 
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obtainable as a degenerate form of Euler’s equation and its adjoint. 


We find that the fundamental solution can be written 


47i U(X, Y; x, y) =f x)*-l(t—X)-*(y/Y) )-wa-o OE (31) 


& 
where C” is now a ‘crossed’ contour passing positively round ¢ = x 
and negatively round ¢= 1. This obeys the ‘characteristic’ con- 
ditions: 
ons = JX, U x y?Xia-X) 
and so satisfies the equation 

8’;U = X(8’—b)U; 
ong= f, U «x (1—z2z)2-1 
and so satisfies the equation 

8U = 2(6—a+1)U; 
atz=JI,y=—f, U = (1—X)-1. 


4. Extension to more independent variables 
The foregoing arguments extend without difficulty to hyper- 
geometric partial differential equations in more than two independent 
variables. Let us begin with the extension of Euler’s equation and 
its adjoint 
88'8"V = ayz(8—m)(8’—n)(8"—p)V 
88'8"U = xyz(8+-m+1)(8’+n+1)(8"+p+1)U 


where 5” = 20/éz. The appropriate characteristic conditions are 


(32) 


these :* 
ona= JZ, 8'8"U = Xy2(8'+n+1)(6"+p+1U, ete.; 
along y= Y, z = Z, 8U = 2YZ(8+m-+1), etc.; 
atz=J,y=Y,2z=2Z, U =(l—iZYZ)"* 

A fundamental solution satisfying these conditions is 

py sree 

ey} \e n! p\(l—ayz)™+1 
This can be written as a repeated contour integral 


(271)? du | dt a Bh de te , (34) 
(l—atu)™+(t—y)"+1(u—z)P + 


* 3, § 6, 302 (top of the page). 


U4, 7,2; 2, 9,2) = ( 
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where C, is a contour in the ¢-plane including ¢ = y and excluding 
t = (xu)-1, and C, is a contour in the w-plane including wu = z and 
excluding wu = a-. 

For the general hypergeometric equation in three independent 


variables f(3)F(8’)4(8")V = xyzg(8)G(8"(3")V (35) 
and its adjoint 

f(—8)F(—8’)¢(—8")U = ayzg(—b—1)G(—8’—1)p(—8"—1)U, (36) 
we write down the fundamental solution 


(onipey = [ [FEY eeatpe tu yt n(ul HOHE Lele Z))}, 
G CG (37) 
where &(x), (x) are associated pairs of solutions of the ordinary 
differential equations 
f(—8)U = xg(—8—-1)U, —f(8)\V = xg(9)V, (38) 
and so forth, and the contours C,, C, are those already defined. 
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J ON THE CONNECTIVITY OF SPACES OF 
POSITIVE CURVATURE* 


By' J. L. SYNGE (Toronto) 
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WE are concerned with a Riemannian space Vy with positive-definite 
line-element, complete and analytic in the sense of Myers.+ It is of 
positive curvature, that is to say, the Riemannian curvature asso- 
ciated with every 2-element (or, equivalently, the Gaussian curvature 
of every 2-space formed by geodesics emanating from a 2-element) 
is greater than some positive number. 

If we take any closed curve C in Vy, the process of parallel pro- 
pagation around C gives a transformation of the vector space at any 
point P on C. By well-known properties of parallel propagation, this 
transformation is orthogonal: thus, if we denote the vector spaces 
before and after by S, S’ and the transformation by 7’, we have 

sg’ = TS, (1) 
where the determinant det 7 is +1. If we repeat the process, 
obtaining the vector space S”, we have 

B= TH = TS, (2) 
where det 7 = (det 7’)? = +1. 

We shall call a curve positively closed if det T = +1, and negatively 
closed if det T = —1. Obviously every curve taken twice over is 
positively closed. The positive or negative character of the closure 
is not altered by deformation of the curve. 

If every closed curve in Vy is positively closed, then Vy is orientable, 
and conversely. 

We shall prove the following theorem: 

THEOREM I. Let Vy be a complete analytic Riemannian space of 
even dimensionality N, positive-definite line-element, and positive 
curvature. Then 

(i) every positively-closed curve in Vy can be continuously contracted 
to a point; 


* Read at the International Congress of Mathematicians, Oslo, 1936. 
+ 8S. B. Myers, Duke Math. Journal, 1 (1935), 40. 
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(ii) every closed curve in Vy taken twice over can be continuously 
contracted to a point; 

(iii) if Vy is orientable, it is simply connected.* 

Parts (ii) and (iii) are obvious deductions from (i): hence it will 
be sufficient to establish (i). This we shall do by proving that, if 
(i) is not true, a contradiction is implied. 

Let C be a positively-closed curve which cannot be contracted to 
a point. (C may be a curve C, taken several times over.) Let L be 
the lower bound of the lengths of all curves into which C can be de- 
formed. We know that Vy is closed;} hence the limit will be attained, 
and there will exist a positively-closed curve G of length Z, such that 
C can be continuously deformed into G. Obviously, G must be a 
geodesic, since otherwise its length could be further reduced, below 
the lower bound L. 

Parallel propagation round G gives an orthogonal transformation 
T of the N-dimensional vector space at a point on G; since G is 
positively closed, det 77 = +1. Let S be the (N—1)-dimensional 
space of vectors orthogonal to G. Since vectors remain orthogonal 
to a geodesic under parallel propagation, this orthogonal vector space 
undergoes an orthogonal transformation 


S’ = TS. (3) 
Since the tangent vector to the geodesic is invariant for the trans- 
formation 7', we have 


det 7 = det T = +1. (4) 


(If G were negatively closed, we should have det 7’ = det 7 = —1.) 

Now an orthogonal transformation in a vector space of M dimen- 
sions leaves a line invariant provided that M is odd: if the deter- 
minant of the transformation is positive, there is an invariant unit 
vector (z’' = p'): if the determinant is negative there is a reversed 
vector (u’? = —yp‘). 

If, as we shall now suppose, the dimensionality N of Vy is even, 
then (VN—1) is odd. Then 7 is an orthogonal transformation of 
positive determinant in a vector space of odd dimensionality. Hence, 
if G is a positively-closed geodesic in a space of even dimensionality, 


* (iii) for N = 2 was proved by Hadamard, but not from the intrinsic point 
of view: J. Hadamard, J. de Math. 3 (1897), 352. 
+ Cf. Myers, loc. cit. 42. 
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there exists at any point P on G a unit vector p' which is unchanged 
as a result of parallel propagation round G. (If G were negatively 
closed, there would be a unit vector reversed as a result of parallel 
propagation round G.) 

Let us now generate a Vy by drawing geodesics g through all the 
points of G in the directions indicated by parallel propagation of 
this ‘invariant’ unit vector yp’, and let us consider the singly-infinite 
family of curves obtained by giving to all the points of G a constant 
displacement v along g. Let us take in Vy a system of Gaussian 
coordinates (wu, v), wu being equal to the arc on G and running between 
the same terminal values (0, Z) on all the curves v = constant, which 
are, of course, closed curves on account of the fact that p* is un- 
changed as a result of parallel propagation round G. 

Comparing the length of G with the lengths of neighbouring curves 
v = constant, we know* that the first variation is zero since the 
displacement is normal to G, and that the second variation is 


BL = }6v7(7T—I1,+ 1,3), i 
where 
T = [4 cos A, 4) J" *, 


L 
I= Ey cos O(E, 7) du, 





0 (6) 
L 
i= {7*s8in?O(€, 7) — Ry jy €' gE" y} du; 
0 J 
here f= a t= = (7) 
ou dv 


&', 7 are absolute derivatives with respect to uw; 7' an absolute 
derivative with respect to v; é, 7, 7 the magnitudes of these vectors; 
6 the angles between vectors indicated; and R&;,,, the curvature 
tensor. Since v is the are on g, and u the are on G, then 7’, €' are 
unit vectors, and since g are geodesics, 7‘ = 0. Also, since n° = p' 
and this vector is propagated parallel to itself along G, we have 
7’ = 0. Hence L 
SL = —}5e* | K ds, (8) 
0 


* J. L. Synge, Proc. London Math. Soc. 25 (1925), 247-64. There is a mis- 
print in equation (6.21); for O(€, n), read O(€, 7). 
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the integral being taken round G, and K being the Riemannian 
curvature of Vy for the 2-element of V,.* Since, by hypothesis, K is 
positive, 5° < 0, and hence the length of G can be reduced by 
variation along g. But this contradicts the assumption that L is the 
lower bound of all curves into which C can be deformed. Hence part 
(i) of the theorem is proved, and, as remarked above, (ii) and (iii) are 
immediate consequences. 

It is unfortunate that the restrictions to positive-definite line- 
element and to even dimensionality deprive the results of interest 
as far as the physical problem of space-time is concerned. For space- 
time has an indefinite line-element, and a space-like section of 
space-time, although it has a positive-definite line-element, is of odd 
dimensionality. It may be of interest to add the following variational 
result, valid for spaces of odd dimensionality, although it does not 
lead to a simple topological consequence. 

Let @ be a (positively or negatively) closed geodesic in a Vy of 
positive-definite line-element and positive curvature: N may be odd 
or even. G may indicate a geodesic G, taken more than once over. 
Starting from a point P on G, parallel propagation round G gives an 
orthogonal transformation of vectors normal to G, a unit vector p‘ 
transforming into pu’. It will in general be impossible to choose p‘ so 
that »’' = p'. But to any choice of y‘ at a point P, there corresponds 
an angle ¢, given by 

cos$ = ay; pp’, (9) 
a,; being the fundamental tensor. If in the flat (N—1)-dimensional 
vector space of vectors normal to G at P we construct ~ the 
(N—2)-dimensional unit sphere with centre at P, then ¢ is the 
length of a geodesic (great-circle) on { joining the points where = 
is cut by the vectors p‘ and p”. 

The angle ¢ will have a lower bound « for all choices of 4‘ normal 
to Gat P. We may call « the ‘twist’ of the geodesic. (If N is even 
and G is positively closed, a = 0.) Let (u')» be the unit vector at P 
which gives this minimum value to ¢, and let it become (u"')p as 
a result of parallel propagation round G. Let us now propagate (,.’')p 
round G in the same sense as before, returning to P with (u”*)p. 


* Since it is proved in the paper referred to above that K is also the Gauss- 
ian curvature of V,, we might have obtained (8) as a variational result 
intrinsic to V,. 
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Now let 
p’ = unit vector field along G, the result of parallel pro- 
pagation of (‘)p, 
pn’? = unit vector field along G, the result of parallel pro- 
pagation of (u”")p. 


(10) 


Then along Gwe have = a,;x‘p") = cosa. (11) 
Let us now define round G a unit vector v‘ by the equation 
vi = [p'sin(sa/L)+p’"‘sin(a—sa/L)|/sin «, (12) 
where s is the arc of G measured from P and L the total length of G. 
Then v‘ is a unit vector which equals (’'), for s = 0 and again (u’')p 
for s = L, because pu‘ arrives at P with the value (y’*),: thus v‘ returns 
unchanged to P. We have then, taking the absolute derivative, 
v' = ofp 'cos(sa/L)—p"'cos(a—sa/L)|/(L sin «), 
p= o2/L2, (13) 
Let us now draw geodesics g through the points of G in the direc- 
tions indicated by v‘ and consider the curves into which G may be 
deformed by taking equal lengths v along g. On account of the way 
in which v‘ has been defined these are closed curves. Passing from G 
to a neighbouring curve v = constant, the first variation vanishes, and 
in the second variation (5) we have, as before, 7’ = I, = 0. We have 


also 2—p—=o2/L?, 0,5) = O(€,v) = 4x, (14) 


and hence the second variation is 
L 


82L = 48? [ (@L-*—K)d (15) 
0 


Hence we have this result: 

THEOREM II. Let G be a closed geodesic of length L in a Riemannian 
space with positive-definite line-element and positive curvature, greater 
than a constant Ky. If the twist « of G is less than LK}, then G can 
be deformed into a shorter closed curve. 

Since a is necessarily not greater than 7, it follows that G can be 
deformed into a shorter curve if 

L > 7K. (16) 
This fact is, of course, well known in connexion with conjugate 
points.* 

* Cf. J. L. Synge, Duke Math. Journal, 1 (1935), 528, for a treatment by 
means of the second variation. 
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